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Abstract 

We consider a non-local spin valve in a Van der Pauw cross geometry with four ferromagnetic 
electrodes. Two antiparallel ferromagnets are used as (charge) source and drain while the de- 
tector circuit involves measuring the voltage between two collinear ferromagnets with parallel or 
antiparallel magnetizations. We find a potentially large increase of the non-local spin voltage. The 
setup displays several additional interesting properties: (i) infinite GMR for the non-local resistance 
(if a symmetry requirement for the device is met); (ii) ON-OFF switch effect, when the injector 
electrodes are parallel instead of antiparallel; (iii) insensitivity to offset voltages. The device can 
additionally be used as a Direct Spin Hall Effect probe and as a reprogrammable magneto-logic 
gate implementing basic operations (NOR, NAND, inverter, AND, OR, etc). 
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I. INTRODUCTION. 



Pure spin manipulation is an important topic in spintronics due to possible applications 
for programmable logic and memory. Non-local spin valves^ - - provide an example of pure 
spin current generation. In the latter a pure spin current generated at a ferromagnetic - 
paramagnetic interface reaches a ferromagnetic probe in the absence of charge current; the 
spin accumulation created in the probe vicinity can then be detected as a charge voltage 
by virtue of Johnson-Silsbee charge-spin coupling 1 ' 7,8 . Such pure spin currents have already 
proved useful to switch magnetizations 9,10 . This is opening promises for further applications 
in logic, sensing and memory devices but to that end it is desirable to increase the signals: 
the spin voltages are typically in the \xV range while mV would be more suitable to ensure 
sufficient SNR (signal to noise ratio). 

While spin valves are in everyday use in hard-drives, in order to reach or even go beyond 
the 1 Tbit/inch 2 density in hard-drives^"—, novel spin valves are required with a low RA 
resistance times area product (typically RA < 0.1 Qfim 2 )^-~— while sustaining a mA current 
and mV voltage; this is beyond TMR (tunneling magnetoresistance) sensors capability since 
they are too resistive. Metallic spin valves are therefore more appropriate. However in order 
for them to have a suitable SNR it is also necessary that the read-heads function with a 
large enough contrast 12 : ARA > 5 mQfim 2 . Metallic spin valves with larger GMR ratio 
AR/R are therefore required. 

At first sight non-local metallic spin valves are not obvious candidates for larger MR 
(magnetoresistance) ratios: they indeed underperform when compared to their local coun- 
terparts; using the same materials and dimensions a local CPP spin valve is expected to 
have a larger AR A since spin confinement is better— —. 

The goal of this paper is to discuss a non-local spintronics device with potentially: 

• enhanced spin voltage in the mV range for currents ~ mA (so that the non-local 
resistance variation AR n i is in the Ohm range) with realistic density currents j < 
10 s A/cm 2 addressing the needs of industry. 

• enhanced non-local GMR ratio: AR n i/R n i (up to 100% for the pessimistic ratio; or 
up to infinity for the optimistic ratio), helping quite generally for better SNR and 
perhaps making them suitable candidates as sensors or read-heads for hard-drive areal 
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densities larger than 1 Tbit/inch 2 . 

Regarding the enhancement of the non-local signal, spin valves with tunnel junctions have 
been reported with non-local resistance in the Q range but due to a polarization decreasing 
rapidly when the current is larger than ~ fiA , the spin voltage remains small in the usual 
fiV range^2. However much progress has been reported recently in pure metallic lateral 
spin valves (~ 10 fiV for nanopillars 10 ) or lateral valves with very low resistance tunnel 
junctions (using a thin nm MgO layer reaching in the latter case the 100 fiV range 

with RA ~ 0.2 Qfim 2 so that already ARA is of the order of a few 1 mQfim 2 . 

We propose to go even further in the improvement by relying on two ideas: (i) use two 
injectors instead of a single one, which should at face value double the signal; (ii) enhance 
the spin confinement by making good use of tunnel barriers, thin enough to stay close to 
the metallic regime but resistive enough to hinder spin leakage. The idea of minimizing the 
spin relaxation volume has been expressed in particular in^ 1 ^ and explains the large signals 
seen in spin valves using carbon nanotubes^. 

Our basic setup applies these ideas by using four collinear ferromagnetic terminals. In 
a standard lateral spin valve^ 1 ^ the charge current flows from a ferromagnetic electrode 
to a paramagnetic drain; the injector electrode is connected by a lateral wire to another 
ferromagnetic electrode used as a detector. In our setup we propose to replace the param- 
agnetic drain by a ferromagnet antiparallel to the terminal acting as current source; the two 
antiparallel electrodes are connected by a paramagnetic metal with thin tunnel barriers in 
order to better confine spin. The two antiparallel ferromagnets act as spin sources although 
in terms of charge one is a source and the other a drain: this effectively doubles the spin 
accumulation in the lateral wire while the tunnel junctions make sure spin is confined. 

We further change the standard detection setup by using a ferromagnetic counter- 
electrode instead of a paramagnetic one. The advantage of using two ferromagnets is 
evidenced when the two detector electrodes are placed symmetrically with respect to the 
injectors (source and drain): provided they are otherwise identical terminals this implies 
that when their magnetizations are parallel, their voltage difference should be identically 
zero by symmetry. This is how we reach an infinite non-local GMR ratio. 

We will also address the issue of voltage offsets plaguing non-local setups^^^ 1 ^: while 
voltages generated by spin accumulation are clearly observed, some additional voltages of 
various origins are also usually seen. These offset (or baseline) voltages have been credited to 
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charge current inhomogeneities2£r— (which impact the calculations done for non-local setups 
since they usually assume one dimensional drift-diffusion equations^i^), or to heating 
(notably Joule and Peltier heating^"—). They may or may not be a nuisance but at any 
rate they prevent observation of pure non-local voltages. The device we discuss in this paper 
can be made insensitive to these offset voltages when the two detector electrodes are identical 
and symmetric since the offsets will cancel out when the voltage difference is measured. This 
is an additional advantage of our device. 

The geometry of our device is that of a Van der Pauw cross as in the Jedema and 
coll. seminal experiments 4 . A close device within a pure lateral geometry will be discussed 
elsewhere^. 

In Section [II] we introduce the Van der Pauw geometry with four ferromagnetic terminals 
and give general expressions for the non-local voltage. The basic functionalities of the device 
are discussed, and notably we will show that the device can perform logic operations (notably 
as a NOR or NAND gate), be reprogrammed to perform other functions (AND, XOR and 
inverter gates), displays a potentially interesting ON-OFF switch effect, and when used as 
a standard 1-bit read-head shows an infinite GMR for the non-local resistance. Use as a 
Direct Spin Hall Effect probe will also be discussed. 

The next section II III studies in detail the impact of the transparency of interfaces and of 
the number of ferromagnetic terminals (two or three out of four) on spin confinement, re- 
sulting in small or large non-local signals. The signals expected are systematically compared 
to those in the standard lateral geometry. 

The last section IIVI discusses the main setup with four ferromagnetic symmetric terminals 
since it displays the previously mentioned properties of (i) immunity to offset voltages and 
(ii) infinite GMR ratio for the non-local resistance. Issues pertaining to the use as a sensor 
are briefly touched upon. 

The bulk of calculations are relegated to the Appendices. Appendix [B] revisits the bipolar 
spin switch calculations by including spin leakage in the measuring electrodes. 
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II. VAN DER PAUW SETUP. 



A. Geometry and notations. 

1. Geometry. 

We consider in this section a four-terminal device in a Van der Pauw geometry (see Fig. 
[1]). The terminals are ferromagnets Fl — F4 positioned as in the Figured! we allow for 
arms of unequal lengths. In sections 111111111 Bl we will allow some of these terminals to be 
paramagnetic through a suitable choice of parameters. 

One-dimensional assumption. We will assume that width and thickness of all arms 
are much smaller than their length. Experimentally, current inhomogeneities due to depar- 
tures from strict one-dimensional flow can arise; however the basic functionalities of our 
device are for the most part independent of that assumption although quantitative predic- 
tions may accordingly lose accuracy. 

Injector electrodes. Fl injects a charge current which is collected in terminal F2. We 
will designate them collectively as injector electrodes; when the need to differentiate them 
shows up, we will say that Fl is the source or injector electrode while F2 is the drain or 
collector electrode. 

Detector electrodes. The detection sub-setup consists in terminals F3 and F4 hooked 
to a voltmeter (or a potentiometer or an ammeter). The latter will measure the non-local 
voltage as a function of the magnetization orientations of each terminal. 

Orientation. We define points 0(x = 0; z = 0) the origin and center of the Van der 
Pauw cross, A(x = 0; z = Li), B(x — 0; z' — L 2 ) , C(x = L 3 ; z — 0) and D(x — 0; z' — L 4 ) 
where each arm has been for later convenience oriented away from O (axis Ox, Ox', Oz and 
Oz 1 ) following Jedema and coll.-. 

The four paramagnetic arms are: / — IV (resp. OA, OB, OC, OD). 



The charge current flowing through Fl — I — II — F2 is I c and flows from top to bottom 
(is therefore negative relative to arm /, but positive, relative to arm //). 
The spin accumulation is defined as: 

(1) 
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Figure 1. Van der Pauw cross with four ferromagnetic terminals. The arm lengths are respectively 
Li, L2, L3 and L4. The arrows represent the magnetization direction. When the device is used as a 
basic spin- valve, only one terminal can switch its magnetization (here F3). The arms are oriented 
away from origin O. 

(where for later convenience we have divided by the electron charge e). 

The spin currents are oriented away from origin O (therefore are counted positive on a 
given arm when flowing away from O). 



2. Spin parameters. 

Arms parameters. The central cross is a normal metal. Its parameters are its spin 
resistance = p* n In/A^ where In is the spin diffusion length, An is the cross-section and 
p* N is the resistivity. We define the lengths of each arm relative to the spin diffusion length 
as (for i — 1 — 4) : 

k = r- ( 2 ) 

Ferromagnet parameters. For each ferromagnetic terminal Fl — F4 (i — 1 — 4), one 
defines the conductivity polarization P Fi (= in Valet-Fert notation 31 ), spin resistance 
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Rf,% = P*f i^F,i/ 'Af,% where lp,i is the spin diffusion length, Af,% is the cross-section and p* Fi = 
(Pit+Pid / 4 • 

(NB: note on terminology; we will call throughout the paper 'spin resistance' the char- 
acteristic resistance found as the product of the resistivity times the spin diffusion length 
divided by the cross-section ). 

Interface parameters. At the interface between the ferromagnets and the paramagnetic 
arms we assume there is a spin dependent interface resistance so that one can define for each 
interface Fi — N (i = 1 — 4) a conductance polarization P cl (= ji in Valet-Fert notation), a 
spin resistance R C i = (Ri^ + Rii) /4. For simplicity we will neglect all spin flips at interfaces 
so that spin relaxation occurs solely in the bulk of the device. 



3. Spin resistance mismatch. 

We define spin resistance mismatch parameters at F - N interfaces as: 

X = (3) 

Three limits can be singled out: 

X < 1: this corresponds to the limit of a transparent junction, which as we will see later 
in detail (Appendix IA2el) is very leaky in terms of spin: this favors large spin currents at 
the cost of reduced spin accumulations in the central paramagnet. 

X > 1: spin confining or tunneling regime, for which spin accumulation increases but 
spin current decreases (in magnitude) (see Appendix IA 2 eh . For X ~ 10, for which the 
contact resistance is moderate (about 10 Q) we will say that we are in the weak tunneling 
limit. This is the most interesting limit in terms of applications to all-metallic read- heads 
or sensors. For R c = 10 2 — 10 4 Q which are usual values in tunnel junctions, X ~ 10 2 — 10 4 
which we will qualify as strong tunneling limit. Although the strong tunneling regime can 
be described by our equations, we will focus primarily in the discussions on the transparent 
and weak tunneling regime where resistances are in the metallic range which interests us for 
sensor applications. 

X = 1: spin impedance matching. The naming for this border situation will be justified 
below in the discussion on effective spin resistance (section III A 5|) . 



4- Effective polarizations. 



The following definition will also prove useful. We define for each terminal (i = 1 — 4) an 
effective polarization as: 

PefU = ^ (4) 

where: 

PRi = (P Fi R Fi + PciRd) I Rn (5) 

and: 

5f = ^f^exp/,±^lexp-/, (6) 
The effective polarization can be rewritten as: 

p _ PFiRFi + PciRci _ 

eff,i ~ R N sinh k + (R Fi + R ci ) cosh k ' ^ ^ 

clearly, |P e //,i| < 1- 

Upon magnetization reversal of the electrode, the effective polarization is an odd function: 

Peff,i > ~ P eff,i- 

In the limit of short arm length 1: 

PpiRpi + PdRci , c s 

— >■ — p — \p) 

KFi + rCci 

which is a weighted average of the electrode bulk and interface polarizations. 

When U — > oo the effective polarization vanishes exponentially which translates the 
complete spin relaxation in the arm: 

PpiRpi ~\~ P c iRci / \ 

Pe f^ = 2 ff +(R I P S eXP ~^ ( 9 ) 

The effective polarization therefore varies between and its maximum value PF ^'^j' c '* R " 
which is bounded from above by sup (Pfi, P C i)- 

5. Effective spin resistances. 

We also define an effective spin resistance for each arm of length U (i = 1 — A) as: 



8 



ReffAki Xj) = Rn "jrr (10) 

_ R Xj cosh k + sinh k _ 
N Xi sinh U + cosh li ' 

R e ff y i(Xi) is an increasing function of spin resistance mismatch Xj. 

It is shown in Appendix IA 2 bl that the spin current l s ,i{0) at the cross center O on arm i 
and the spin accumulation there are related by A/i(0) = —R e ff^ I St i(0) (for arms III — IV; 
for arms I — II a, more general relation taking into account the spin injection at Fl and F2 
holds). This is an analog of Ohm's law for spin which explains our identification of R e ff,i as 
a spin resistance. (Note that the analogy is not complete: the relation for spin is a local one 
(expressed here at point O), while Ohm's law holds for a voltage difference and is therefore 
non-local. This results of course from the non-conservation of spin current.) 

We also define a total effective spin resistance for the device: 

Reff = ^ (12) 

2^=1-4 R efftl 

The previous expression admits obvious generalization to an arbitrary number n > 4 of 
arms. 

As can be seen from its definition, the total effective resistance R e ff is related to the arms 
spin resistances R e ff,i by the analog of a parallel resistance addition law. We will show later 
(section [II B 2j) that the spin voltage is proportional to the total effective spin resistance R e ff 



so that large effective resistances are desirable; this will also effectively demonstrate for our 
geometry the parallel addition law for spin resistances. (For a discussion of spin resistance 
addition law at nodes we refer the reader to^.) Since the total effective spin resistance is 
the sum of four resistances in parallel, whenever one is much smaller than the others, it 
will short the other arms: spin leakage will be stronger so that spin accumulation will be 
reduced. 

The length dependence of the effective resistance for various values of the spin resistance 
mismatch X is shown in Fig. |2]|3] . At large distance the effective resistance converges 
exponentially fast to i? n the spin resistance of the paramagnetic arm: 

Reff,i — > Rn (13) 
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which reflects the fact that the spin relaxation is dominated by the paramagnet bulk. In 
the limit U — > 




Rfi + Rci 



(14) 



(which is sensible since the paramagnet is then too short for spin relaxation to occur). 
When li <C 1, the effective resistance remains close to Rpi + R C i if the distance I obeys: 



When X ^> 1, the condition becomes ! < 1/X which reflects a steeper exponential decrease. 
The effective spin resistance for a given arm is therefore comprised between Rn and 



(or the reverse inequality if Rn > Rf% + Rci)- 

As a rule the effective spin resistance will be larger for large interface or ferromagnet spin 
resistance (R c and Rf)', since the ferromagnet spin resistance Rp is in general much smaller 
than the paramagnet spin resistance Rn due to short spin diffusion lengths, large interface 
resistances R c are required to achieve large effective resistances R e ff,i- 

For X = 1, one observes that the effective spin resistance does not depend any more on the 
arm length and is equal to Rn- One then has (on arms III or IV) A/i(0) = —Rn l s ,i{0) 
which is the relation one would get from an infinite arm. Everything happens as if the 
interface had been washed away: this is the reason why we qualified the case X = 1 as 
corresponding to spin impedance matching in III A3! 

In the transparent regime {X < 1), the effective spin resistance is larger at large distance, 
which is an interesting feature for the design of large non-local circuits. This advantage is 
circumvented by the exponential decrease of the effective polarization (the non-local resis- 
tance will be shown to be proportional to both in III B 21) so that in terms of large signals 
the transparent regime is not interesting, neither in the short-distance nor the large- distance 
limit. 



I \X 2 - l| < X. 



RF,i + Rc,i'- 



Rn < Reff,i < RFi + Rci 



(15) 
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Effective resistance Retf(x, I) 




Distance I 

'o 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 2. Effective resistance for one terminal (normalized to Rn) as a function of arms length I 
(relative to paramagnet spin relaxation length) for various impedance mismatches X = 1 — 10 in 
the (weak) tunneling regime. 

B. General expression of the spin voltage in Van der Pauw geometry. 

1. Spin accumulation at cross center. 

Solving the one- dimensional drift-diffusion equations (see Appendix |A] for details of the 
calculations) leads to the following results. 

The spin accumulation in the center of the cross is found as: 

A/i(0) = (P effA - P e//>2 ) R eff I c (16) 
We can define a total effective polarization at injector electrodes Fl and F2: 

Pinjector = Peff,l ~ Peff,2- (17) 

The minus sign in front of P e ff.2 comes from the fact that F2 is a charge drain: in terms 
of spin accumulation it therefore acts contrariwise to electrode Fl. This also shows that to 
ensure maximum signal it is better to have opposite orientations for Fl and F2 (antiparallel 
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Effective resistance 

1 



X=1 




Distance I 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 3. Effective resistance for one terminal (normalized to Rn) as a function of arms length I 
(relative to paramagnet spin relaxation length) for various impedance mismatches X = 0.1 — 1 in 
the transparent regime (X < 1). 

injector electrodes). This is easy to understand: when injector electrodes are parallel, F2 
acts as a spin sink for the spins injected by Fl; therefore the spin accumulation should 
decrease. But when F2 is antiparallel to Fl, spin leakage is frustrated; although F2 is a 
charge drain, it acts as an additional spin source. 

It is noteworthy that the spin accumulation does not depend on the magnetization orien- 
tations of detector electrodes F3 and FA (although it does depend on its parameters through 
-R e //). This is sensible: spin injection is ensured by Fl and F2 not F3 and FA. 

The general structure of the spin accumulation in terms of Fl and F2 relative orientation 
can be understood simply. Suppose we flip all spins of the setup. Then: A/i — > —A/i 
since we have exchanged spin up and spin down electrons. This implies that the spin 
accumulation must be an odd function of electrode polarizations. This is easily checked 
on Eq. f fl~6l) : when both Fl and F2 are flipped, their effective polarizations get reversed 
P e ff — > —P e ff and therefore the spin accumulation is reversed. Another way to reach the 
same conclusion is to notice that the physics of the setup should be invariant when both 
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current and magnetizations are reversed. Flipping all spins or reversing the charge current 
should lead to the same spin accumulation, namely a reversed one. 

An important consequence of the structure of Eq. (I16p is that the spin accumulation at 
the cross center vanishes when Fl and F2 are identical and parallel so that P e ff,i = Peff,2- 
This follows clearly from symmetry: when Fl and F2 are symmetric with respect to the 
cross center and have identical parameters, the spin accumulation on the line F1 — F2 should 
be antisymmetric and the spin current should be symmetric when Fl and F2 are parallel 
(the opposite when they are antiparallel). This can be understood by reversing the current: 
if one reverses the current, on the one hand, the spin accumulation at the cross center 
should not change since Fl and F2 are identical and parallel (by watching Fig. [TJ after a tt 
rotation); on the other hand, reversing the current must reverse the spin accumulation since 
reversing the current is equivalent to reversing all spins. The only way out is for the spin 
accumulation at the cross center to vanish. This result does not depend on the assumption 
of one-dimensional flow and follows directly from symmetry. 

We will see later that this enables an ON-OFF switch function onto the device. 

2. Spin voltage and non-local resistance. 

Spin voltage. The spin voltage (or non-local voltage) is the voltage drop between 
terminals F3 and F4: 



The behaviour under magnetization reversal is easy to understand: when the magnetiza- 
tions of both detector electrodes are switched, their coupling to the spin accumulation gets 
reversed so that the non-local spin voltage should change sign. The spin voltage is therefore 
odd under magnetization switching of both detector electrodes. 
Inserting the expression of the spin accumulation in Eq. ( fl6l ): 




(18) 



Straightforward calculations (see Appendix IA2 dl) lead to: 



V nl = -A/i(0) (P e// , 3 - P effA ) . 



(19) 




(20) 



This expression factors out neatly in three contributions: 
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(i) the geometry dependent effective spin resistance Reff = (^=1-4 Reffi) > 

(ii) a total effective polarization for injector and collector electrodes (Fl and F2) 



injector e//,l 



-P. 



eff,2, 



(21) 



(iii) and a total effective polarization for the two detector electrodes (F3 and FA): 



F detector — Peff,3 ~ Peff,4- 



(22) 



Non-local resistance. The non-local resistance (sometimes called a transresistance) 
is defined as the ratio of the non-local voltage to the charge current flowing through the 
injector electrode Fl to the collector electrode F2 : 



To achieve a large signal it is therefore necessary to have a large total effective spin resistance 
for the cross R e ff and to have on the one hand antiparallel source and drain terminals, on 
the other hand antiparallel detector electrodes. 

R e ff is largest when all arms effective spin resistances are also large, which is the case if 
spin resistance mismatches are in the tunneling regime according to the discussion in III A 51 
and if the arms length is short enough. Two situations may arise: 

i) one or several spin resistance mismatches are in the transparent regime (X < 1). Then 
R e ff ~ Rn or smaller. No enhancement of the spin voltage is to be expected when compared 
with the usual lateral setup. We will say that we have an open geometry^ which leaks spin 
(larger spin currents but smaller spin accumulations). 

ii) all spin resistance mismatches are in the spin confining regime {X > 1). Then R e ff ~ 
R c in the limit of short length for the arms (and assuming for simplicity mismatches roughly 
equal Xi ~ X ~ R c /Rn- In such a geometry which will be qualified as closed^ the signal is 
therefore potentially much larger than in an open geometry. 

Let us compare the non-local resistance to local resistances in the same device, 
i) Firstly a local resistance can be measured between source and drain (Fl and F2); as 
shown in Appendix IA 3 at 




(23) 



Reff [Peff,l - Peff,2f + Rq {Pl-z} 



(24) 



14 



where 

i=l-2 

,p i [PFjRciRFi - PcjRli - ZPFiPciRFiRci] / ot -\ 
ripi + -Kci 

(zi are the locations of probes in the electrodes, see Appendix IA 3 ap . In general Ri oca i,i2 ~ Rq 
since the GMR effect is a few percents; the non-local resistance which is commensurate with 
Riocai,u — Ro is therefore much smaller than Ri OC ai,i2- 

It is more meaningful to compare the variations upon magnetization switching AR = 
Rap — Rp] for the non-local signal, we have set in the following Fl and F2 antiparallel 
while switching F4 magnetization: 

ARnl = (Peff,l + Peff?) PeffA ^ 
^Rlocal,12 2Peff,l Peff,2 

(where for Ri oca i, either terminal Fl or F2 have been switched). In the case of identical 
source and drain electrodes (Fl and F2) this reduces to: 

ARnl = PeffA . ^27) 
ARl cal,l2 PeffA 

if the ratio is larger than unity (-P e //,4 > PeffA) there is an amplification of non-local MR 
variation versus local MR (or the converse if P e ffA > PeffA)- 

ii) If we then compare to the local resistance found when current flows from F3 to F4, 
one gets instead: 

(28) 



AR nt _ (PeffA + Peff,2) 



&-Rlocal,M 2P e //,3 

One can again get an MR amplification (or reduction if P e ff,3 is sufficiently large). 
But if all electrodes are identical, all these MR ratios are then equal to unity 

ARnl 



(29) 



^Rlocal 

which means the non-local measurement performs equally well as local measurements in 
terms of raw resistance variation, with one proviso: non-local resistances have smaller base- 
lines. Indeed for R n i, the baseline or smallest signal is: 



effAl 



\Rnl ,min 

I = R eff |(|P e //,l| - |-P e // )2 |) (|-Pe//,3| ■ 

while for Ri oca iA2- 

\Rlocal,min\ = Ro {Pl-2} + Reff (|-P e //,l| — \Peff,2\) 2 
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which is larger on account of the Ro term. 

iii) It is worthwhile to compare to CPP GMR for a spin valve with similar dimensions: 



m N (p c Rc+p f Rf) 2 



2R N (R F +R c )coshl+[(R F +R c ) 2 +R 2 N 'j sinhi 

for identical and infinite ferromagnetic layers separated by the same distance L = I 1^. In 
terms of X = (Rp + R c ) /R^ this can be recast as: 

a D &Rn {RcRc + PfRf) 2 
IXricpp — n o 

(x + iy expi- (x - iy exp —4 

so that for a non-local device with identical terminals: 

ARnl 1 (30) 



^Rcpp 2 

The non-local signal is smaller by a factor 2 which is easy to understand: in the cross 
geometry the spin relaxation volume is doubled when compared with a spin valve with a 
paramagnetic layer of identical length because of the side arms. 

A systematic comparison of the non-local resistance in the cross geometry with the stan- 
dard lateral geometry is left to sections IIII AlUVl 



3. Non-local charge current. 

The voltage between the detector electrodes F3 and FA actually acts as an electromotive 
force (emf) of magnetic origin; when F3 and FA are shorted, a charge current therefore 
appears. This non-local charge current is induced by the spin accumulation generated by 
the remote source and drain Fl and F2. Note that since the spin voltage is an electromotive 
force, it can be measured either through a (nano-) voltmeter or with a potentiometer: in the 
latter case, the advantage is that there is no current at all during the measurement. 

If however it proves advantageous that the signal be a current (for chaining the non-local 
device to a bipolar transistor for instance rather than a MOSFET), the non-local current is 
found (see Appendix IA 3 bp as: 



Inl = h (31) 

where i?/ ca/,34 is the local resistance measured when one drives a current between F3 and 
FA. This can be rewritten as: 

Reft (Peffl - P ef f,2) (PeffS ~ I'.fJ,) , 



Reff [Peff,3 - PeffA} 2 + Ro {-^3-4} 
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(32) 



where: 

^0 {iV4} = P*Fi i 1 - P Fi) Z i +PNU + R* 

i=3-4 

_|_ [PFjRgRFi — PcjRli ~ ZPFiPciRFiRci] 
RFi + Rci 

(it can be checked that Rq > 0); Z{ (i = 3 — 4) are the locations on F3 and F4 of the wires 
which short them together (in the following we have chosen to focus on R n i rather than I n {). 



C. Main properties of the device. 

The device functionalities depend on the terminals for which the magnetization has been 
fixed (prevented from switching); additional symmetry requirements can also add properties. 
We first discuss functionalities pertaining to sensing or data storage. 

As explained in the introduction the non-local resistance and spin voltage are contami- 
nated by offsets, the origin of which is still under debate (current inhomogeneities; thermal 
origin). This may adversely affect the measured signals and in the following we will take 
care to indicate the potential impact of these offsets. 

1-bit reading. This is the basic functionality of the spin valve as a sensor. Suppose the 
orientations of all terminals but one (say F3) are pinned (for definiteness we assume: Fl :j, 
F2 :4- and FA :|. ) so that F3 acts as a sensing electrode. We make no special assump- 
tions on the terminals (later on some conditions will be imposed for further functionalities). 
One recovers a spin-valve behaviour (see Fig. HJ) with two distinct values of the non-local 
resistance which uniquely determine the orientation of terminal F3: 

Rp = Reff {Peff,l + Peff,2) ( P eff,3 ~ P eff,i) 

R AP = R eff (P effA + P e//>2 ) (P e// , 3 + P effA ) 

(the P/AP index refer to F3 magnetization orientation relative to F4). The spin valve can 
therefore be used to read a 1-bit information. 

In the standard non-local spin valve and if we neglect voltage offsets, the spin voltage 
changes sign when one terminal is flipped so that for antiparallel and parallel alignment 
Rap = —Rp (this is recovered here in the limit f e //,4 = 0, when FA is a paramagnet); this 
is not the case here although depending on the relative values of P e ff,3 and P e ff,A there can 
still be a change of sign (Fig. 0]-b, when P e ff,3 > Peff,i)- 
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Figure 4. Non-local resistance variation when the orientation of a single terminal (here F3) is 
switched by an external magnetic field (here #3). All other terminals are pinned. Injector Fl 
and collector F2 are antiparallel with orientation as shown on the graph, (a) case P e ff3 < Peff,A'- 
standard spin valve effect (1-bit reading); (b) case P e ff,3 > Peff,4' s P m valve effect with spin 
voltage change of sign.). 

For the Van der Pauw cross under the most general conditions the difference with the 
standard non-local setup is therefore minor; it remains to see if larger signals can be achieved. 
This will be the topic of sections IIIIHIVI where we will show that tunnel contacts at the four 
terminals can greatly enhance the non-local resistance. 

Furthermore this functionality is clearly affected by offset voltages which will shift the 
signals and change the GMR ratios (adversely if the offset is positive). We now discuss a 
simple way to circumvent these offsets. 



Offset free 1-bit reading with infinite non-local GMR. 

The offset issue is easily fixed if we assume F3 and FA are identical electrodes placed 
symmetrically with respect to the rest of the setup. Due to symmetry offset voltages are 
neutralized (assuming offset voltages are spin independent) since they will shift both voltages 
V3 and V4 in the same manner so that the spin voltage V n \ = V3 — V4 is free from offsets. 

Additionally: Rp = R e ff (P e ff,i + Peff,%) (-Peff,3~Peff,4) = in parallel alignment 
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since -P e //,3 — Peff,A while R n i 7^ in antiparallel alignment (we have assumed that Fl and 
F2 are antiparallel and pinned). 

One has achieved an infinite GMR for the non-local resistance since the ratio GMR = 
(Rap — Rp) I ' Rp — > oo (using the optimistic ratio; the pessimistic ratio would be 100 %). 
There is an intrinsic contrast which is protected by symmetry from the voltage offsets. 

Such a maximized MR ratio is clearly helpful for SNR in terms of Johnson noise (or shot 
noise in the strong tunneling limit) since the latter scales as AR/\/R~; indeed if we compare 
with a CPP spin valve with a AR/R =1 — 10% GMR ratio, this would imply at identical 
AR (AR = AR local = AR nl ) an increase of SNR by 10 - 20 dB. (Indeed: SNR for non-local 
device would be oc 1/V AR so that S N R non _[ oca i / SNRcpp oc a/ Riocai / AR) . The difference 
is quite significant given that under operation one expects in general at least 30 dB SNR. 

In terms of geometry requirements, note that the symmetry between detector electrodes 
is required only on the scale of a few spin relaxation lengths lp (on the ferromagnet side) 
since the spin accumulation is washed at larger distance. 

The property is also clearly independent of the precise geometric arrangement, does not 
depend on the assumption of one-dimensional flow and will be valid for other geometries 
than the cross studied in this paper, provided the two detectors are arranged symmetrically 
with respect to the injectors. Experimentally this is very useful since this gives a lot of 
flexibility in terms of design. 

3-bit reading or storage. For that function one terminal is pinned, while the other 
three are free and play the role of input signals (the non-local resistance measured between 
F3 and FA is as previously the output signal). Eq. ([23]) shows the non-local resistance can 
assume 8 different values when the electrodes orientation are changed and the maximum 
value for \R n i\ is reached when on the one hand Fl and F2 are antiparallel, and on the 
other hand F3 and F4 are also antiparallel. The orientations of the three non-pinned 
terminals are uniquely determined by the 8 distinct values of the spin voltage. This implies 
that the device encodes 3 bits in principle (3 bit spin valve). 

Fig. [5] shows the eight outputs signals as a function of the input variables (for illustrative 
purposes F3 is varied by an external field H3 on the graphs, showing the output variation 
with the change of one input). 

The property survives offset voltages which come as an additive contribution to the 
voltage. However resolution may be adversely affected by offsets. Note that although there 
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Figure 5. Non-local resistance variation as a function of terminals magnetization orientation showing 
3-bit sensor or storage. For illustrative purposes F3 magnetization is varied in each figure (a-d) 
(through an applied field H3) to show the output signal variation with the change of one input, 
(a) Injector Fl and collector F2 are antiparallel. (b) They are switched to parallel: the signal 
collapses (to some extent), (c) Signal is reversed when injector and collector are both switched 
from antiparallel alignment in (a), (d) Injector and Collector are parallel (but with directions 
opposite to (b) ). 



are four ferromagnets and a priori 4 bits could be stored, the spin voltage can assume only 
8 values, not 2 4 = 16. This is because when all spins are reversed, the spin voltage is 
unchanged in Eq. ( 1231) . 



Offset-free 2-bit reading. By the same token applied previously to detectors for the 
1-bit reading it is possible to convert the 3-bit reading function into a protected 2-bit reading 
function: suppose that detector terminals F3 and F4 are symmetric (namely, have identical 
parameters and are placed symmetrically with respect to the device) but that Fl and F2 
differ (P e //,i 7^ Peff,2)- Let us pin terminals 3 and 4 in antiparallel orientation (F3 :j and 
F4 while only Fl and F2 are allowed to switch magnetizations. The non-local spin 
voltage then takes 4 distinct values depending on the orientations of terminals 1 and 2: 
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R n i = R e ff (P e ff,3 + Peff,4) (±-P e //,i ± Peff^)- These values uniquely determine the 2-bit 
state of Fl and F2. The main advantage when compared with the 3-bit function discussed 
previously is that since terminals F3 and FA are symmetric, there can be no offset voltages 
since they automatically cancel out when measuring the voltage drop between F3 and F4 
(if offset voltages are spin-independent; if it is not the case, they will add up to the spin 
voltage). Note that if Fl and F2 are pinned while F3 and FA are not pinned but identical 
(to avoid offsets), the spin voltage only assumes three values since R P = 0; to store 2 bits 
without offsets it is therefore necessary to pin F3 and FA, not the other way around. 

ON-OFF switch for the 1-bit read-out or storage function. When source Fl 
and drain F2 are identical ferromagnets (same distance from origin, same conductivities, 
polarizations, interface resistances, spin diffusion length), P e ff,i = Peff,2- This then implies 
that when Fl and F2 are parallel, the spin accumulation at the cross center vanishes 

Afi(0) = (34) 

so that whatever the orientation of F3 and F4, Rp = = Rap- The spin voltage has been 
killed and we have disabled the read-head or 1-bit storage. The property is clearly interesting 
in terms of logic if the device is chained to another device for instance a MOSFET whose 
gate is controlled by the spin voltage (after suitable amplification). 

That property should survive offset voltages since symmetry protects it. 

But one might wonder if the spin voltage measured at F3 and F4 will still vanish if we 
take into account departures from strict one-dimensional charge and spin flow. It is clear 
indeed that even if we take into account the 3-dimensional nature of the device but remain 
in a quasi-one-dimensional approximation, the spin accumulation on the side arms (zones 
III — IV) will remain small (it may be non zero at edges) and will never diffuse as far 
as F3 and F4 provided the width of each arm is much smaller than its length (w <C I): 
indeed when Fl and F2 are parallel, the spin accumulation is an odd function of z in the 
direction Fl — F2; therefore, on the side arms (FiF — IV), there will be some spilling of spin 
accumulation with opposite signs on opposite edges, close to the origin O. But if w <C I, 
spin diffusion will mix these opposite spin accumulations which will cancel out when one 
reaches the detector terminals (F3 and F4). The spin voltage measured at F3 and F4 
should therefore still vanish. 

Direct Spin Hall Effect probe. The geometry lends itself easily to probing the Spin 
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Hall Effectp^ 1 ^. Imagine Fl and F2 are normal electrodes and that the magnetizations of 
F3 and FA are perpendicular to the plane of the cross. On the arm AB (Fl — I — O — 1 1 — F2 
or zOz') a spin accumulation may appear due to Hall coupling an = o- xx /a xy on the width 
w (for —w/2 < x < w/2): 

Afi(x, z) = aH x. (35) 

When we reach the cross center, there can be a spilling of charge current lines to the side 
arms, but we will neglect this effect by assuming that the arms width are much smaller than 
their lengths (w <C Li for i = 1 — 4). When we move on the side arms (/// — IV), the spin 
accumulation decreases in magnitude. This decrease can be estimated using drift-diffusion 
equations as: 

A/i m (C) = ^ Af,(x = z = 0) (36) 

S 3 2 



see Appendix IA 2 dl Eq. (IA48I) ] with a similar expression for arm IV: 



Aii iv (D) = ^ Afi(x = ~,z = 0) (37) 
Oj 2 



A///7j(C) can be re-expressed as: 

5+ L x + L 2 2 



A/^m(C) = — - - (38) 



and: 



Afi IV (D) = -— -. (39) 

dj L\ + L 2 2 



The spin voltage is therefore (see Appendix IA2dl) : 



w 

V„, = - (Pe//,s + P. tM ) Aa(x = y) (40) 
= -(P^ + P^^-l (41) 

Note that for identical and antiparallel electrodes F3 and FA, the signal therefore vanishes 
since P e ff,3 = —P e ff,4 and that for maximal non-local signal, parallel magnetizations are 
required. 

Let us go back to our initial setup with four ferromagnetic electrodes. If magnetizations 
are perpendicular to the plane of the Van der Pauw cross, the spin Hall effect will add up to 
the non-local spin voltage arising from spin injection. But if we choose a symmetric setup 
with identical parallel electrodes F3 and FA, offsets cancel out in V n i — V 3 — V4 and Vp — 

22 



as explained previously. Therefore the only remaining signal is that of the Direct Spin Hall 
Effect. This provides an interesting all-electrical alternative to the observation of Direct 
Spin Hall Effect, which initially was observed optically in semiconductors^ 1 ^ although later 
on some electrical detection schemes have been used in metallic systems to investigate both 
Direct and Inverse Spin Hall Effector— . 

D. Implementation of magneto-logic gates. 

Many proposals exist in the literature for the use of magnetoelectronics circuits as logic 
gates 45 ; this has prompted a lot of activity in the field of semiconductor spintronics since 
integration to existing processes would be optimal 4 ^. Our device is metallic and therefore 
not the best candidate as a spin transistor since there is no amplification: this renders 
the chaining of gates more delicate for instance (unless one uses hybrid designs combining 
pure magnetoelectronics devices with conventional transistors). However, the device still 
possesses obvious capability as a programmable magneto-logic gate as we now demonstrate. 

Let us associate bit with down 4- magnetization and bit 1 with up magnetiza- 
tion. Fl and F2 are assumed to be identical; F3 and F4 are also identical (same arm 
length, same polarizations, spin diffusion length, etc); only the magnetization orientations 
are allowed to differ. The non-local resistance can then assume only three values: and 
±i?o = i4i? e jj P e ff,i Peff,3- For logic operations we will consider two conventions: 

(i) associate bit to zero resistance, and bit 1 to R n i = +Rq. 

(ii) or (opposite convention) associate bit 1 to zero resistance, and bit to R n i = +R®. 
Note that in what follows we have discarded all configurations for which R n i = —Ro to 

avoid ambiguities in the bit association and keep only those for which there are only two 
possible outputs: and Rq. 

For the first convention assigning the bit content of Rni, the expression for R n i [see Eq. 
f l2"3"j) ] can therefore be rewritten in terms of bits as the Boolean equation: 



(Fl - F2) (F3 - FA) = R nl 



(42) 



where to simplify notation we have conflated R n i and its bit content. 
With the opposite convention, 



(Fl - F2) (F3 - F4) = R, 



(43) 
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Note that in order to have well defined HIGH and LOW states the associated voltages (or 
non-local resistances) must be sufficiently different which is not ensured in the presence of 
offset voltages. That's why we choose symmetric terminals to get rid of these offsets. 

We first try to reproduce basic binary Boolean functions. We need therefore to pin two 
terminals; the other two terminals will represent the variables treated by the device in the 
following manner (for instance): imagine each terminal is screened to prevent any magnetic 
field applied to a given electrode to influence any other electrode; we then apply external 
magnetic fields to either of the two non-pinned terminals but not to the other twos, whose 
magnetizations are therefore fixed during the whole operation. 

NOR gate. We choose convention of Eq. ( l4"2l) for bit coding. We pin F2 and F4 (F2 :| 
and F4 :f) so that F2 = 1 and F4 = 1 in bit terms. To have a non-zero resistance among 
the four possible bit configurations of Fl and F3, only Fl = and F3 = are allowed. 
The non-local resistance is then R n i = +Ro which we associate with HIGH state or bit 1: 
this means that Fl F3 = R n i; this can also be recovered through algebra by using Eq. ( )42l) 
which in our case is: 

(Fl - 1) (F3 - 1) = R nl (44) 

so that: 

FlF3 = R nl . (45) 

This is precisely a NOR gate since Fl F3 = Fl + F3 by De Morgan theorem. This is 
a very important property since NOR has functional completeness: any Boolean function 
can be implemented by using a combination of NOR gates (only NAND possesses the same 
property). 

OR gate. In the same configuration as for the NOR gate, if LOW state (bit 0) is now 
associated with R n i = +Ro and HIGH state (bit 1) to R n i = [convention of Eq. ([4*3]) ]. one 
obviously gets an OR gate. Of course since the association of HIGH and LOW has been 
reversed with the previous case, the two settings are incompatible since they correspond to 
opposite bit assignment. 

AND gate. We still pin F2 and F4 (F2 :| and F4 :|) so that F2 = and F4 = in 
bit terms. Using Eq. ( T42l) this implies Fl F3 = Rni which is an AND operation (using the 
first convention for R n i bit content). This gate is compatible with the NOR gate but not 
the OR gate which is produced with a different convention for bit coding. 
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NAND gate. By changing the bit coding of R n i, the AND gate turns into a NAND gate 
with the same configuration for F2 and FA. As mentioned previously, the achievement of a 
NAND is quit noteworthy since it has functional completeness in Boolean algebra. 

AB or A B (A not implied by B) gate. We now pin F2 and F3 and choose F2 = 1 
and F3 = in bit terms. Then Eq. (H2l) (first convention) turns into Fl FA = R n i. 

A + B or A -<= B gate. Using the same configuration for F2 and F3 as previous gate 
but exchanging the conventions for the bit content for R n i leads to an inverted gate: indeed 
Fl FA = R n i implies Fl + FA = R n i by De Morgan theorem. 

AB or A =£> B gate. We pin F2 and F3 and choose F2 = and F3 = 1 in bit terms. 
Then Eq. ( H2l) (first convention) turns into Fl FA = R n i. 

A + B or A =>- B gate. Choosing the second convention of Eq. ( H3l) with F2 = and 
F3 = 1 leads to Fl FA = R n \ which is equivalent to Fl + FA = R n i by De Morgan theorem. 

FALSE gate. We pin Fl and F2 in parallel configuration. Then R n i = (in terms of 
resistance not bit) whatever the configuration of F3 and FA. If we adopt the first convention, 
the zero resistance translates into bit 0. This is therefore a FALSE gate. 

TRUE gate. In the same terminals configuration as previous gates, if we adopt the 
opposite convention for the bit content of R n i, one then gets a TRUE function. 

The six other possible binary operations (out of sixteen) can not be built as easily using 
a single cross device; this leaves out the NOR and XNOR gates from the list of the basic 
gates in use in electronic logic (AND, OR, NOR, NAND, XOR, XNOR, buffer, inverter). A 
way to achieve them would be to chain gates since all Boolean functions can be recovered 
using only a NOR (or NAND) gate. 

We now turn to the two other basic gates, buffer and inverter. 

Buffer. For that function only one terminal is not pinned. We choose to pin Fl = 1, 
F2 = and FA = 0. Using Eq. ( H2l) this implies F3 = R n i. This can also be realized by 
pinning Fl = 0, F2 = 1 and FA = 1. Using Eq. (S3D this implies F3 = R^. 

Inverter. When pinning Fl = 1, F2 = and FA = 0, using Eq. (03]) implies F3 = R n i. 
This can also be realized with Eq. ( H2l) by pinning Fl = 0, F2 = 1 and FA = 1 which leads 
toF3 = F n ,. 

The NOR, AND, A B, A =fr B , FALSE, buffer and inverter gates use the same 
bit convention for R n i (HIGH state for R n i = +Rq, LOW state for R n i = 0) and can be 
implemented together in the same circuit; since they differ only by the assignment of pinned 
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terminals, this means that one can easily turn a gate into another one (reprogrammable logic) 
through local application of an external magnetic field or through spin transfer torque. 

In summary, the non-local cross with four ferromagnetic terminals is a versatile spintronics 
device which can be used either as a standard spin valve: when some symmetry requirements 
are met, it displays the important property of an infinite GMR for the non-local resistance; 
it can be immunized against offset voltages observed in many non-local setups. Finally, it 
can be used as a magneto- logic gate as well as probe the Direct Spin Hall Effect. 

III. VAN DER PAUW CROSS WITH TWO OR THREE FERROMAGNETS. 

Which principles should guide us in order to achieve large signals in non-local setups? 
As stressed time and again in the literature^^ 3 -, they are quite simple: larger spin accu- 
mulations in the paramagnet can be generated if (i) the paramagnet volume in which spins 
can relax is small (in comparison with l^) and if (ii) spin back-flow to the ferromagnets is 
hindered by large enough interface resistances. 

One can classify geometries as open or closed according to the (non-) fulfillment of these 
prescriptions 23 : it has been argued in the latter reference that at small enough volume and 
for large enough tunnel barriers the spin accumulation can be greatly enhanced because it 
then scales with the (large) interface resistances R c . But whenever the geometry is open, the 
spin accumulation scales with the much smaller spin resistance of the paramagnetic channel 
i?7V which in practice is often in the Ohm range. Although prescriptions (i) and (ii) have 
been stressed repeatedly, spin leakage in the current or voltage probes is often overlooked 
although they may alter significantly the signal (an example is provided in Appendix IBl). 

The goal of this section is to examine the transition from open to closed geometries as 
a strategy to enhance the non-local signal used for instance to read or sense 1 bit, and to 
study its interaction with the doubling of spin injector electrodes. The signal generated in 
our cross geometry will be systematically compared with the spin voltage in the standard 
lateral device. 

The following general features will be a guide: a large signal requires a large total effective 
spin resistance for the device R e ff- Since the latter can be interpreted as resulting from the 
addition of four spin resistances R e ff,i (i = 1 — 4) in parallel corresponding to the four arms 
of the cross, one will expect a small signal whenever one or several of these spin resistances 
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are significantly smaller than the others since they will short the other arms. If all effective 
spin resistances R e ff,i are commensurate, in order to have a large R e ff,i, it is then necessary 
that R ci ^> Rj\j simultaneously for all four terminals as explained in section Hi B 21 

(As an aside remark we wish to bring to the reader's attention that the expressions derived 
in the literature^ 1 ^ for the bipolar spin switch transistor of Johnson remarkably show the 
enhancements characteristic of closed geometries over open geometries with a typical scaling 
with interface resistance R c at small enough distance between the ferromagnetic terminals. 
As discussed in Appendix [B] this is due to neglecting spin leakage to the current drain; when 
this leakage is taken into account, the signal is actually found to scale with R N as in open 
geometries.) 

A. Van der Pauw cross with two ferromagnets. 

1. Open geometry with two ferromagnets. 

For the sake of comparison we first consider the case when the detector electrodes are 
paramagnets as in the cross arms / — IV. When Fl and F3 are identical ferromagnets at 
the same distance of origin (li = l 3 ) while F2 and F4 are identical paramagnets with spin 
resistance Rn, the spin resistance mismatches at each terminal are 



X 1 — X 3 — X; X 2 



X 4 = 1 



(46) 



while: 



PR 2 = PRi 



PRi = PR 3 







PR 



(47) 
(48) 



and defining 5p for arms / and 




exp li ± 



(X-l) 



exp 



(49) 
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exp Z 2 / 4 



(50) 



27 



so that: 

2 



R, 



ni 



R N (PR 



a l( r 3 R N (PP)' 



— — (51) 
[(1 + X) 2 exp/ + (X 2 -l)] 

where I = 2L\/l^ the total length separating the ferromagnetic terminals (in units of the 

spin diffusion length in the paramagnet) and where we have defined a± = ±1, 03 = ±1 to 

index the majority spin directions of electrodes Fl and P3 (relative to an absolute axis). 

In the standard geometry, since one has a spin valve it is customary to quote the resistance 

variation when one ferromagnet magnetization is switched; this is twice the maximum value 

Rp\ 



2R N (PR\ 

5R ri i n = Rp — Rap = t n r- (52) 

' [(1 + X) 2 exp/ + (X2-l)] 1 ; 

This generalizes the expression found in the literature for the non-local resistance variation: 
more precisely the result quoted by Jedema and coll.- corresponds to the case of vanishing 
interface resistance, so that X = Rf/Rn {— 1/M using Jedema and coll. notations 4 ). 

Although the geometry is open (in the sense that spin current can leak easily since 
terminals P2/P4 do not hinder its flow (X 2 = X 4 = 1) ), it is interesting to observe that the 
larger the resistance mismatch X at the ferromagnets, the larger the signal (for instance, 
if we set P = 1 in Eq. (JH2D one gets that 5R ntfi ex X 2 / [(1 + X) 2 expZ + (X 2 - 1)] which 
is an increasing function of X). This means that large resistance mismatches are already 
beneficial and increase the spin accumulation although there is some spin leakage. 

The largest value is at short distance when the denominator of Eq. (152]) is 2X (X + 1). 
In the tunneling regime X 1, 5R n ifi cx R^. In the opposite limit Z < 1, the signal will 
be even smaller since it scales as R c + R F R N : 

6R nlfi (l — )> 0) > inf (P c , P F f (R c + R F ) 

and 

8R nl;0 {l —►())< sup (P c , P F f {R c + R F ) . 

Such a scaling which is at most ~ R^ or even below (<C Rn) is characteristic of open 
geometries. 
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Let us compare to the non-local resistance variation for the lateral geometry with iden- 
tical parameters (same distance I between injector and detector, same spin resistance mis- 
match X at ferromagnets, same cross-section for the paramagnet channel connecting the 
ferromagnets) 6 : 

4-R P~R 

SRnl,lateral = r 2 " 1 IT- (53) 

[{2X + 1) exp / — exp —/J 
We plot on Fig. |6] the ratio of the cross signal versus the one in the lateral geometry: 

fY n SRnifi [(2X + l) 2 exp/-exp-/] 

m{X, I) = — ; — = — p 7) —. (54) 

" ttfil, lateral Z \(1+Xy exp/ + (X 2 -l) 

In the short distance limit: 

m(X, I — > 0) — ► 1 

which is expected since the cross and the lateral geometries are then identical. At large 
distance, the signal is larger in the cross geometry whenever X > l/\/z~ (which includes 
the tunneling regime at the ferromagnets and also part of the transparent regime). The 
impact of a large value of X is moderate (about 10 % between X = 10 — 100) in stark 
contrast to what we will observe in a closed geometry. The general behaviour is easy to 
understand: in the tunneling regime spin current leakage is less pronounced in the cross 
geometry since the paramagnetic drain is further away (while the detector ferromagnet and 
paramagnetic counter-electrodes are at the same distance). In the transparent regime, by a 
similar reasoning one gets the opposite. 



2. Spin confining geometry with two ferromagnets. 

It has been argued by Jaffres and coll.— that spin confinement tends to increase non-local 
signals (resistance or voltage) since spin accumulation is stronger whenever spin leaking is 
hindered by large tunnel barriers. We can study this by considering that electrodes F2 and 
F4 are normal paramagnets but that there is a resistance mismatch 

Y = (55) 

it AT 
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Figure 6. Non-local resistance in cross setup compared with lateral setup. I is the distance in units 
of spin diffusion length between injector and detector in either setup. In tunneling regime (X > 1) 
and in part of the transparent regime, the signal is larger in the cross geometry. 

where Rn,o and R c ^ are the spin resistance of F2 and FA and the interface resistance between 
them and the central cross. One can therefore interpolate between an open geometry (the 
standard geometry, i.e. Y = 1) and a closed geometry (y>l and X ^> 1). We will now 
study the transition from one regime to the other. 

Let us assume that terminals Fl and F3 are identical ferromagnets at identical distance 
li = l 3 = l of the cross center O, that F2 and FA are identical paramagnets also at the 
same distance l 2 = U = Iq from O {ABCD is then a square); the total distance from the 
injector Fl to the detector F3 is I = 2l . 

We also set: 



Xi — X 3 — X; X 2 — X 4 — Y 



(56) 



and: 



PRx = PR 3 = PR, 

PR 2 = PR 4 = 0. 



(57) 
(58) 
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There are several obvious ways to create this spin resistance mismatch Y: one is to 
deposit a tunnel barrier between the terminals F2 and F4 and the central cross (R c ,o 7^ 0); 
another is to use the same paramagnet for both F2, F4 and the central cross but have a 
different cross-section (R Ci o = but Rn,q = P*n^n/^n,o 7^ Rn = Pn^nI-^n)- 

One ends up with: 

a,a 3 R N (PR) 2 

Rnl ~ 



°N 



2 F 0p + & 

where 



and 



so that: 



A 2F 

The signal will be largest at short distance for which: 



, 2 -.. XY 



and: 



inf (P c , P F y R n — — < 6R nl (l — )• 0) 
A + Y 



XY 

8R nl (l^0) <sup(P c , P F f R N 



X + Y 



so that 8R n i(l — > 0) scales as R^ x+y 

Therefore in the strong tunneling regime (X ^> 1 and Y ^> 1), 5R n i scales as 

inf (X, Y) R N > R N 



(59) 



X± _ r± r± _ r± 
°F ~ w l/3' W JV ~~ °2/4- 

[The definitions for 8f are given in Eq. ([6]).] 
Let us define 

A 2F (X, Y,l)= 2 {(Y + 1) exp l + (Y- 1) exp -/q]- 1 

{(Y+l) exp Z [(X + l) 2 exp 2l + (X 2 - 1) 
- (Y - 1) exp -l [(X - l) 2 exp -2l + (X 2 - 1)] } (60) 



A (X, /) = [{2X + l) 2 exp 2/ - exp(-Z)] (61) 



4i?7v (Pi?) 2 
= ^ ^ (62) 



yielding much larger signals than in the lateral geometry. But whenever either of X or Y is 
in the transparent regime, 5R n i will scale as the smaller of the two and so will be at most at 
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the scale of Rjy as predicted for open geometries (see III B"2j) . In order to check the impact 
of increasing the spin resistance mismatches at the various terminals to the signal in the 
lateral geometry we consider a non-local resistance ratio m(X, Y, Z) per: 

m(X, Y, I) = r j,'> (63a) 

Orinl, later al\^i L ) 



Ap(X, I) 
A 2F (X, Y, I)' 

The relative increase of the signal is then m(X, Y, I) — 1. At / = 0: 



(63b) 



Y (X + l) 

m(X, Y,l = 0)^ ^ + y } . (64) 



At large distance: 



m(X, Y, I) — > (2X + l) \ (65) 
V ; 2 (X + l) 2 V ' 

which is larger than 1 whenever x > 1/ \/2 and more generally belongs to the interval [0.5; 2]. 

The large distance behaviour is of course less interesting since the non-local signal is weaker, 

so we will discuss in detail only the short distance behaviour. 

At short distances, large values are achieved whenever both parameters X and Y are large 

(spin confining regime, closed geometry). For X ^> Y ^> 1, m(X, Y, I — 0) ~ Y; while for 

Y > X > 1, m(X, Y, I = 0) ~ X+l; finally for X ~ Y » 1, m(X, Y, Z = 0) ~ (X + 1) /2. 
This third situation is probably the easiest to achieve: indeed if we want to gain at least an 
order of magnitude the other cases would imply that one of the mismatches is at least two 
order of magnitudes larger while when X ~ Y 3> 1 they both have the same magnitude. 
For a realistic value X ~ 10 in the weak tunneling regime (where the contact resistance 
is usually metallic) this yields an enhancement by a factor up to 5 (see Section [IV] for a 
discussion of realistic parameters X ~ 10). 

We plot in the next figures (Fig. [TH!]) the length dependence for various values of X and 

Y in the weak tunneling regime (keeping moderate values (X, Y) < 10) covering the three 
situations described in the previous paragraph (larger X , larger Y or equal X = Y). 
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Non-local resistance ratio m(X, Y, I) 




Figure 7. Ratio of non-local resistance variation in cross versus lateral setup as a function of 
distance I between injector Fl and detector F3 (in units of spin diffusion length Ijy) for spin 
resistance mismatches Y = 10 (at terminals F2/FA) and X = 1 — 10 (at terminals F1/F3). Fl 
and F3 are identical ferromagnetic metals while F2 and FA are paramagnets in this section with a 
spin resistance mismatch with the central paramagnet. 



Still in the weak tunneling regime we also plot on Fig. [10] the non-local resistance ratio 
for fixed distance I = 0.1 — 0.3 as a function of spin resistance mismatch X (assuming 
Y = X). Values in the range I = 0.1 — 0.3 are quite reasonable experimentally (for instance 
for Cu at ambient temperature Zjv = 300 nm while a distance L = 100 nm is within reach 
of lithography). At such distances the enhancement can still be several hundred of percents 
as can be seen in Fig. [10] unless X < 1 (transparent regime). 



The transparent regime is less interesting in terms of an increase of the signal since the 
geometry is now open. Yet whenever Y is larger than 1, even if X is in the transparent 
regime, one still gets an increase of the signal at short distance (Fig. ITT|) . But if both 
parameters are smaller than unity, the signal gets reduced by several tens of percent when 
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Figure 8. Ratio of non-local resistance variation in cross versus lateral setup as a function of 
distance I between injector Fl and detector F3 (in units of spin diffusion length In) for spin 
resistance mismatches Y = 1 — 10 (at terminals F2/F4) and X = 10 (at terminals F1/F3). Fl 
and F3 are identical ferromagnetic metals while F2 and FA are paramagnets in this section with a 
spin resistance mismatch with the central paramagnet. 

compared against the lateral setup as can be seen on Fig. [T2]fT3l 



To summarize this section devoted to the cross geometry with two ferromagnets (one as 
charge source and spin injector, the other as spin accumulation detector), we have confirmed 
the impact on spin confinement of tunnel barriers even of moderate strength (X = 10 might 
correspond to R c ~ 10 Q since Rn is usually in the Ohm range). In the transparent limit, 
the device under-performs when compared to the standard lateral spin valve and should be 
avoided. 
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Non-local resistance ratio m(X, Y, I) 




Figure 9. Ratio of non-local resistance variation in cross versus lateral setup as a function of 
distance I between injector Fl and detector F3 (in units of spin diffusion length ijy) f° r s pi n 
resistance mismatches Y = 1 — 10 (at terminals F2/FA) and X = Y (at terminals F1/F3). Fl 
and F3 are identical ferromagnetic metals while F2 and FA are paramagnets in this section.. 

B. Van der Pauw cross with three ferromagnets. 



We now use two ferromagnets as spin injectors. When the (charge) source and drain 
are in antiparallel orientation the signal is enhanced because both ferromagnets acts as spin 
sources. But when they are parallel, one is a spin source and the other a spin sink so that 
spin accumulation is reduced. If the electrodes are identical ferromagnets, then the signal 
will be doubled when compared with the case where there is only one spin injector electrode. 
This is exactly what we found since R n i oc P e ff,i ~ Peff,2 = 2P e // 5 i if electrodes are identical 
and antiparallel. 

We assume that Fl — Fl — F3 are identical electrodes (Fl and F2 antiparallel while F3 
can switch from one orientation to the other) and FA is a paramagnet. The spin resistance 
mismatches are set as: 



Xi — X2 — X% — X; X4 — Y. 
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(66) 



Non-local resistance ratio m(X, Y=X, I) 




Spin resistance mismatch X 
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Figure 10. Non-local resistance ratio in cross versus lateral setup at fixed distances I = 0.1 — 0.3 
between injector Fl and detector F3 (in units of spin diffusion length Zjv) for spin resistance 
mismatches X = Y. 



We also assume h = h = h = h an d use the same definitions as in previous section: 

PR 1 = PR 2 = PR 3 = PR 
PR 4 = 



(67) 
(68) 



Then: 



R 



2(7ia 3 R N [PR 



(69) 



where o\ = ±1 and 03 = ± refer to the majority spin direction of ferromagnets Fl and F3 
relative to an absolute axis. 

Let us define 

A 3F (X, Y, I) = [2(Y+ 1) exp l + 2 (Y - 1) exp -k}' 1 

x {(V + 1) exp l [2 (X + l) 2 exp 2/ - (X - l) 2 exp -2l + (X 2 - l)] 
- (Y - 1) exp -la [2 (X - if exp -2l - (X + l) 2 exp 2l + (X 2 - l){}0) 
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Figure 11. Ratio of non-local resistance variation in cross versus lateral setup as a function of 
distance I between injector Fl and detector F3 for spin resistance mismatches Y = 10 (at terminals 
F2/F4) in the weak tunneling regime and X = 0.1 — 1 (at terminals F1/F3) in the transparent 
regime. 



Then: 



5R 



ni 



Rp — Rap 



AR N (PR 



A 



3F 



(71) 



where P/AP refer to the direction of F3 relative to Fl. Therefore the non-local resistance 
roughly scales as: 



5R 



C .R/v ^ 2 



ui 



(72) 



AaFpr,, y, /) 

(where C is a constant; we have used upper and lower bounds on ( Pi?) as in section UlIAl 
above). 

The largest signal is found at small distance when: 

XY 



5R 



R 



N 



3Y + X 

We observe a characteristic quadratic dependence of the numerator against a linear one in 
the denominator as a function of spin resistance mismatch: this is what ensures that in the 
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Figure 12. Ratio of non-local resistance variation in cross versus lateral setup as a function of 
distance I between injector Fl and detector F3 for spin resistance mismatches Y = 0.5 (at terminals 
F2/F4) and X = 0.1 — 1 (at terminals F1/F3). All electrodes are in the transparent regime. 

tunneling limit, 5R n i can scale as R c (here as the smaller of R^X ^> Rn or RnY 3> R^. 
But as soon the geometry is open due to one or both spin resistance mismatches in the 
transparent limit, 5R n i scales as R n or below. 

The ratio of the non-local resistance with regards to the standard lateral setup is then: 

A 



m(X, Y } I) 



(73) 



By varying the spin resistance mismatch Y at terminal F4 one can interpolate between 
an open geometry (Y < 1) and a closed one (Y ^> 1 while I>1), 

Let us make some general comments comparing the three ferromagnets cross setup with 
the two ferromagnets cross discussed in the previous section IIIIA2I 

When X — Y, one gets A 2 = 2A 3 so that 

m 3F (X, Y = X, l) = 2m 2F (X, Y = X, I). 

The non-local signal is exactly twice that found for the setup with two ferromagnets. Indeed 
when all terminals have the same spin resistance mismatch, the effective resistances R e ff are 
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Figure 13. Ratio of non-local resistance variation in cross versus lateral setup as a function of 
distance I between injector Fl and detector F3 for spin resistance mismatches Y = 0.1 — 1 (at 
terminals F2/F4) and X = Y (at terminals F1/F3). All electrodes are in the transparent regime. 

identical in both setups since the parameters are identical; but the total effective polarization 
at injector -P e //,i — Peff,2 = 2Peff,i is doubled in the 3 ferromagnet setup (-P e //,2 = in the 
two ferromagnet setup). 

For X > Y the effective resistance (see III B II) for the 3 ferromagnet setup becomes larger 
than that of the 2 ferromagnet setup (see IIII A~2~| this is quite clear since R e ff^(Xi) is an 
increasing function of Xi so that 

3 1 



eff,3F _R eff>1 (X) + R eff)1 (Y) 

r 2 2 

Reff,2F ~ S m + 



> 



Reff,l{X) R efftl (Y)_ 

This means that the non-local signal will be more than doubled with respect to the 2 
ferromagnet setup studied in IIII A 21 This effect does not require a totally closed geometry 
to occur. 

For X < Y, the 3 ferromagnet geometry is by the same arguments less spin confining 
since we have three terminals with a smaller spin resistance against only two in the 2 ferro- 
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magnet setup. However the doubling of effective polarizations remain which may mitigate 
the decrease of effective spin resistance. But as a rule in order to achieve stronger signals 
one should seek the condition X > Y. 

1. Open geometry (Y = 1). 

We first consider the case when the paramagnetic terminal has no tunnel barrier but is 
perfectly matched in terms of spin resistance to the central cross. This is in the terminology 
of Ref.— an open geometry. Nevertheless as in the 2 ferromagnet geometry studied in IIII Al 
we will find that spin accumulation can still be enhanced even though not all the terminals 
are tunnel barriers. 

The spin resistance mismatch Y is set to unity: X4 = Y = 1. 

The non-local resistance ratio to the two ferromagnet cross (with Y = 1) is then: 

m(X, Y = 1,1) =2 [{2X + 1) 2 exp I - exp -l) (74) 

x [2 (X + if exp I -(X - if exp -/ + (X 2 - l)] -1 (75) 
At large distance its limit is: 

m(x, y = i, 0-+ i^irj ■ (76) 

which is always larger than 1 whether X is in the spin confining (X > 1) or the transparent 
regime (X < 1). 

At small distance the ratio tends to: 

m(X,Y = l, Z)^4^±l (77) 

which is larger than 2 and can be as large as 4 in the tunneling regime (provided X > 1) 
and is comprised in the interval [4/3; 2] for X < 1 (transparent regime). 

So although the geometry is open, the use of tunneling junctions does have some spin 
confining effect, which here is reinforced by the use of two injector electrodes resulting in an 
enhancement of the non-local resistance which can be up to 300% increase in the tunneling 
regime. For a distance I = 0.2 and X = 2 — 10, m = 2.45 — 3.30 (145 — 230 % increase). 
As can be seen in the next figure (Fig. [T4"I) the enhancement can be seen at all distances so 
that this geometry is already better than the lateral one in terms of spin confinement. 
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Figure 14. Ratio of non-local resistance in cross geometry (with three ferromagnets) versus lateral 
setup as a function of distance I between injector and detector. Spin resistance mismatch is X = 
1 — 10 at terminals Fl — F2 — F3 (weak tunneling regime). 



In the transparent regime {X < 1), the doubling effect coming from the two spin injectors 
compensate partly the spin leakage due to transparent junctions so that there is always at 
least a 33 % increase at short distance. The enhancement is still present at large distance 
and can then reach up to 125 % as can be seen in FigJT5l 



2. General case. 

We now conjugate the effects of spin confinement by tunnel barriers and the doubling of 
injector terminals. At small distance one gets: 

m(X, Y = 1,1) — y AY X + l (78) 
v , > J X + 3Y y ' 
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■local resistance ratio m(X, Y=1, 1) 




Figure 15. Ratio of non-local resistance in cross geometry (with three ferromagnets) versus lateral 
setup as a function of distance I between injector and detector. Spin resistance mismatch is X = 
0.1 — 1 at terminals Fl — F2 — F3 (transparent regime). There is still an enhancement of the signal. 



•_'.v • r 2 



while at large distances: 

m ( X , Y — 1, I — y oo) — >■ [^j^ ) < 4- (79) 

Large values are achieved whenever both parameters X and Y are large (both in the 
tunneling regime). 

Let us focus first on that tunneling regime. 
There are then three interesting limits: 

(i) X ^> Y ^> 1, which implies m(X, Y, I = 0) ~ 4Y; (note also that the spin confinement 
is more pronounced in that limit than in the 2 ferromagnet setup due to a larger effective 
spin resistance); 

(ii) Y > X > 1 implying m(X, Y, I = 0) ~ f (X + 1); 

(iii) finally for X ~ Y > 1, m(X, Y, I = 0) ~ X + 1 (see Fig. QMS]). 

For a realistic value X = Y ~ 10 (see section HVCj) this yields an enhancement by a factor 
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Figure 16. Ratio of non-local resistance in cross geometry (with three ferromagnets) versus lateral 
setup as a function of distance I between injector and detector. Spin resistance mismatch is X = 10 
at terminals Fl — F2 — F3 and Y = 1 — 10 at paramagnet FA. All electrodes are in the tunneling 
regime. 



up to 11 reaching therefore an order of magnitude. At a finite distance the enhancement 
remains considerable ranging from m = 3.3 — 6.6 at I = 0.2 for X = 10 and Y = 1 — 10 (Fig. 
161). 



If one exchanges X and Y (with large Y — 10 and X = 1 — 10), one gets weaker signals. 
This is normal, since when Y is larger than X, spin confinement is disfavored since there 
are then three terminals with spin resistance mismatch smaller than the fourth (Fig. H7|) . 



When X and Y are about equal and large, X ~ Y ^> 1 the relative increase (~ X + 1) 
for a given mismatch X is smaller than when F^X^>l(m~| (X + 1)); however the 
latter condition might be more inconvenient to realize for a modest gain (for instance, if we 
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Figure 17. Ratio of non-local resistance in cross geometry (with three ferromagnets) versus lateral 
setup as a function of distance I between injector and detector. Spin resistance mismatch is varied 
{X = 1 — 10) at terminals Fl — F2 — F3 and set at Y = 10 at paramagnet FA. All electrodes are 
in the tunneling regime.. 

aim at an order of magnitude increase, X ~ 10 this would imply Y ~ 100 (in the strong 
tunneling limit) yielding m ~ 14 at I = 0; while already for X = Y ~ 10, m ~ 11 at I = 
(see Fig. QSJ. 

In practice the distance between injector and collector can not be reduced arbitrarily. 
Nevertheless the ratio m (X, Y, I) may remain large as can be seen in the next figure (Fig. 
[T9|) which shows the relative increase of signal at distances I — 0.1 — 0.3. For X ~ Y = 10, 
m ~ 5.7 — 8.0 for I = 0.1 — 0.3 (which are quite reachable at low temperature for I = 0.1 or 
even room temperature for / = 0.3). This is a doubling of the signal when compared to the 
values we found for X ~ Y = 10 in the two ferromagnet closed geometry (of section IIII A 21) . 
(As previously mentioned, this stems from the fact that the total effective polarization is 
doubled due to the use of two spin injectors while the effective resistances are identical.) 
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Figure 18. Ratio of non-local resistance in cross geometry (with three ferromagnets) versus lateral 
setup as a function of distance I between injector and detector. Spin resistance mismatch is varied 
X = 1 — 10 at terminals Fl — F2 — F3 but Y = X at paramagnet FA. All electrodes are in the 
tunneling regime.. 



Turning now to the transparent regime, we still observe an enhancement of the signal 
(see Fig j20|) : as before this is entirely due to the doubling of spin injector terminals since 
there is spin leakage at all terminals. For X — Y < 1, the ratio m varies between X + 1 and 
2 for an increase up to 125 % at large distance. 



To summarize this section, we have seen the impact of doubling the number of spin 
injectors conjugated to spin confinement resulting from large spin resistance mismatches. 
Although this geometry with three ferromagnets is quite interesting in terms of the large 
signals which can be achieved, in practice it is better to consider a four ferromagnet setup 
which will have all the qualities of the 3 ferromagnet Van der Pauw cross with some additional 
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Figure 19. Ratio of non-local resistance in cross geometry (with three ferromagnets) versus lateral 
setup at distances I = 0.1 — 0.3 between injector and detector. Spin resistance mismatches are set 
to identical values at all terminals X = Y 

properties: protection from voltage offsets (which may be detrimental to SNR) and infinite 
GMR ratio for the non-local resistance. 



IV. VAN DER PAUW CROSS WITH FOUR FERROMAGNETS. 

We now consider now the main device of this paper, a four ferromagnet Van der Pauw 
setup with the following symmetries: (i) identical injector and collector electrodes (Fl and 
F2), (ii) identical detector electrodes (F3 and F4). Besides an enhancement of the spin 
voltage such a device is protected against voltage offsets and offers an infinite non-local 
GMR ratio. Other general properties are described in section Hi CI notably such a setup will 
also display an ON-OFF switch effect (for the basic 1-bit basic sensing function). 

46 



2.4 
2.2 

2 
1.8 
1.6 
1.4 
1.2 

1 

0.8 
0.6 



Non-local resistance ratio m(X, Y=X, I) 



X=1 



X=0.7 



X=0.5 



X=0.4 



X=0.2 



distance I 



Figure 20. Ratio of non-local resistance in cross geometry (with three ferromagnets) versus lateral 
setup as a function of distance I between injector and detector. Spin resistance mismatches are set 
to identical values at all terminals X = Y in the transparent regime (X = 0.1 — 1). 

A. Non-local resistance and infinite GMR. 

We consider a rhombus geometry for the cross: OA = OB = l\ while OC = OD = Z3. 
The resistance mismatch parameters for the device are set as: 



while: 



Xi — X2 — X; X3 — X4 — Y 



PRi = PR 2 = PR, 



PR* = PRa = PR 



(80) 

(81) 
(82) 



The non-local resistance [Eq. 02 3|) ] becomes: 



Rnl — Reff (fl — 0"2) (°3 — 04) Peff,l Peff,3 



(83) 



where <jj = ±1 refer to the majority spin direction of terminal i relative to an absolute axis. 
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Eventually: 

R Rn (<n - <J_2) (<r 3 ~ ^4) PR , . 

nl [(X + l)(F + l)exp/-(X-l)(F-l)exp-/] 1 ' 

where 

/ = Zi + h- (85) 

Upon magnetization switching of one detector electrode the variation in non-local resistance 
is therefore: 

4R n PR PR 

nl \{X + 1) (Y + 1) exp I - (X - 1) (Y - 1) exp -/] 
B. Amplitude of non-local resistance. 

The signal is largest at small distance where the non-local resistance obeys: 

XY 

inf (P c , P F ) inf (i* P' F ) AR n — — < 5R nl (l — )• 0) (86) 

l A + Y ) 

with a similar upper bound. Again one recovers the characteristic quadratic dependence on 
spin resistance mismatches at the numerator (linear at denominator), which imply that if 
both X and Y are large, SR n i will scale with them, much beyond R^. If one of the spin 
mismatches (say X) is in the transparent regime, it will set the scale, below or at most at 
Rn- 

Let us compare in details with the standard lateral setup; we will set PR = PR and 
X = Y for simplicity (When I/F the same trends against the case X = Y are seen as 
in the three ferromagnet geometry.) 

Then the non-local resistance ratio becomes: 

m (X I)- 6Rnl - [(2^ + 1) 2 exp /-exp-/] 

SR n l,lateral [(X + l) 2 exp / — (X — l) 2 exp — /] 

At small distance, one gets the enhancement: 

m(X, 1)—>X + 1. (87) 

This is the same enhancement as in the 3 ferromagnet setup when X = Y (see IIII B~2j) . 

We first consider the tunneling regime (X > 1). We plot in the following figures (Fig. 
(I21H22P ) the length and resistance mismatch dependence of the ratio m(X, I). For X — 1, 
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Figure 21. Ratio of non-local resistance in cross setup (with four ferromagnets) versus standard 
lateral setup, as a function of distance I between injector and detector. Spin resistance mismatch 
X is in the spin confining regime (X = 1 — 10). 

we observe that m ~ 2 at all distances: in the case of open cross geometry with two 
ferromagnets (see section UlI A lj) already one had m ~ 1; the doubling here stems therefore 
from the two spin injectors. For X = 10, the enhancement remains strong at increasing 
distance: for instance for I — 0.1 — 0.3, m ~ 500 — 770 %. 



For completeness we now turn to the transparent regime (X < 1) although it is less 
interesting than the tunneling regime if one is to achieve large signals. We observe the same 
trend in Fig. [23] as before with an increase in several tens of % and up to 100 % at large 
distance, which stems from the double spin injectors. At shorter distance, the spin current 
leakage is more pronounced. 
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Figure 22. Ratio of non-local resistance in cross setup (with four ferromagnets) versus standard 
lateral setup, at fixed distance I = 0.1 — 0.3 between injector and detector. Spin resistance mismatch 
X > 1 is in the spin confining regime. 

C. Discussion and comparison to experiments. 

Can we expect signals in the mV range for non-local voltages? Already the state-of-art 
has reached ~ 100 — 200 fiV in lateral spin valves^; the peculiarity of these lateral spin 
valves is that they use thin tunnel MgO barriers (1 — 6 nm) between the ferromagnets 
(Permalloy NiFe) and the paramagnetic channel which is silver Ag. Due to the thinness 
of MgO the interface resistances are much smaller than in usual tunnel junctions, with an 
interface resistance times cross section product in the range RjA ~ 0.2 Qfim 2 . While for as 
deposited samples the signals are already higher than usual (~ 10 fiV), after annealing the 
signals can jump by an order of magnitude; this has been accounted on oxygen vacancies 
which increase after annealing. 

Let us see how much we can get in similar conditions but with the cross geometry with 
four ferromagnetic electrodes: we will borrow the parameters found of Fukuma and coll.— 
for lateral spin valves based on Py — MgO — Ag junctions; the size of the junctions is 
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Figure 23. Ratio of non-local resistance in cross setup (with four ferromagnets) versus standard 
lateral setup, as a function of distance I between injector and detector. Spin resistance mismatch 
X is in the transparent regime (X = 0.1 — 1). All terminals are identical in this plot. 

A = 0.022 /im 2 ; the separation between detector and injector is L = 300 nm . 

At 10 K, for a 1 mA current, the largest spin voltage is measured at 112 fiV with the 
following parameters: P c = 0.44 for the interface conductance polarization (noted Pj by 
Fukuma and coll. 21 ), Pp = 0.35 for Py electrodes, lp = 5 nm, = 1100 nm for the spin 
relaxation lengths of Py and silver Ag, Rn = p* N l^/A = 0.89 Q for Ag spin resistance, 
Rf = pplpjA = 0.09 Q {RNiFe = 0.08 Q in Fukuma and coll.— is related to our spin 
resistance Rp by RNiFe — Rf (1 — Pp) ), Rc — 12.13 Q (our interface spin resistance R c 
is related to Ri of Fukuma and coll. by R c = Ri/(l — Pj) with R t A = 0.2152 fi/im 2 and 
A = 0.022 fim 2 ). 

Then the spin impedance mismatch is as large as: 

X = = 13.73 

Rn 

while: 

I = L/l N = 0.27. 
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For these values and using the measured AV n i = 110 fxV, AR n i = 110 mfl, 1 = 1 mA 
this implies: 

m{X, I) = 660 % 

which would imply a spin voltage variation as large as: 

AV ~ 725 fiV 

or: 

AR nl ~ 725 mft 

(since 7 = 1 mA); if we use the symmetric four-terminal cross device, the baseline signal 
vanishes (for parallel orientation of detector electrodes), so that 

R A = AR A = 16 mO/im 2 

which compared with CPP GMR spin valves is large in terms of AR A product (a few 
mQfim 2 ) but small in terms of customary R A values (10 — 100 mQ/xm 2 )^. 

Still at 10 K, the highest voltage reported is 200 fxV for a larger current J = 3.5 mA 
so that the non-local resistance is smaller AR n i = 63 mQ; using our setup we predict one 
would have achieved signals in the range: 

AVni ~ 1-320 mV 
AR nl ~ 415 mVt 
AR A rsj 9.1 mttfim 2 
j ~ 1.6 10 7 'A/ cm 2 . 

The current density is therefore reasonable although the mV mark for the spin voltage is 
reached! The non-local resistance variation is however smaller (Fukuma and coll.— attributes 
the decrease to an enhanced spin relaxation in Ag due to phonons). 

If we put into perspective these predictions, we may remember that the first measurements 
of non-local resistance in bulk metal wires^ by Johnson and Silsbee yielded ~ nQ, were 
thereafter improved in thin films^ to about 10 fiQ; Fukuma and coll. 21 have pushed to 
~ 100 mVt. The largest non-local resistance variation observed is however already in the 
Ohm range in lateral tunnel junctions^ 1 ^, with one major catch: currents must remain 
small (~ ix A) due to spin depolarization at larger currents, which results in a spin voltage 
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still in the customary \lV range. In contrast, in metallic lateral spin valves larger currents 
can be achieved^i^ 1 ^"— . 

Turning now to room temperature, there is few change in spin impedance mismatch X 

X = 13.44 

but due to a much smaller ~ 300 nm 

I = 1. 

One still finds a sizable enhancement: 

m(X, I) = 415 % 

which would lead to the following numbers (using original experimental values of Fukuma 
and coll. 21 AV nl = 51 fxV, AR nl = 51 mfi, I = 1 mi): 

AV n i ~ 210 fjV 
AR nl ~ 210 mfi 
AR A ~ 4.5 mVLfim 2 
j ~ 5 10 6 A/ cm 2 . 

There is still room for improvement at room temperature before reaching the mV mark. 
One obvious way would be to decrease the distance between injector and detector; for 
instance suppose L = 100 nm which is quite reasonable in terms of lithography, then 
/ = 0.33; taking the previous room-temperature value for X (X = 13.44) would yield: 

m(X, I) = 605 % 

so that 

AV n i ~ 310 fiV 
AR ni ~ 310 mtt 
AR A ~ 6.8 mttum 2 
j ~ 5 10 6 A/ cm 2 . 

If we lower to a still achievable 50 nm for which I = 0.17 then: 

m(X, I) = 795 % 
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and: 

AV nl ~ 405 fjV 
AR nl ~ 405 mfi 
Ai? A ~ 8.9 mfi/im 2 
j ~ 5 10 6 A/cm 2 . 

(at such length L = 50 nm m = 1150% at 10-fT so that AV n i ~ 1.250 mV again in the mV 
mark at / = 1 mi). 

Are there any other ways to further improve the signal at room-temperature? We will 
show elsewhere 36 that with an all-lateral structure even stronger enhancements are found 
due to better spin confinement. In the latter example one would get as much as m ~ 1070 % 
at room temperature (or m ~ 1320 % at 10 K). 

Application to sensor and read-heads technology? 

Reaching the 1 Tbit/inch 2 mark up to 10 Tbit/inch 2 areal density requires new recording 
methods to push beyond the super-paramagnetic limit; bit patterning and thermally assisted 
magnetic recording are currently investigated^. But such large areal densities are also very 
challenging for current sensor technologies. The transition from AMR, CIP GMR, TMR 
and now CPP GMR has accompanied dramatic and steady increase of the areal density of 
hard-drives up to about 0.5 — 0.7 Tbit/inch 2 (as of late 2011). Further increase to reach 
the 1 Tbit/inch 2 has motivated the latest transition from TMR read-heads to CPP metallic 
spin valves because TMR is limited to R A ~ 1 Qfim 2 which will not be small enough 
to accommodate the high data rates required for the desired areal densities 12 . In contrast 
CPP GMR read-heads can have R A < 0.1 Qfim 2 which is the minimum required for an 
areal density 1 Tbit/inch 2 . An additional requirement for CPP spin valves is to have larger 
AR A > 5 mQfim 2 instead of the usual< 1 mQfim 2 which are too small to ensure good 
SNR (signal to noise ratio). The source of noises are usually Johnson thermal noise in 
metallic spin valves (both CIP and CPP) or shot noise in TMR valves. But at the shrinking 
sizes relevant for larger areal densities additional relevant sources of noise appear to be STT 
(spin transfer torque) noise, mag-noise (noise due to thermal fluctuations in the free layer) 
and amplifier noise^ 2 -^ 1 ^. 

Non-local devices have been barred for practical use due to small signals and have mostly 
been interesting for fundamental and pedagogical use. As we have shown this may prove less 
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of a concern by using geometries with better spin confinement and with multiple injectors. 
Such non-local devices used as sensors would hold several advantages over their standard 
counterparts while retaining the small RA of CPP metallic spin valves: 

- smaller thickness of the sensing detector (two or three layers since the reference layer is 
located elsewhere and can have a different coercive field by playing on its geometry); 

- diminished sensitivity to STT noise since the largest momentum transfer is between 
source and drain electrodes; 

- no self- field since no current flows through the detectors (Oersted fields are at injectors). 

- In the case of a symmetric cross geometry, the possibility to achieve 100 % or maximum 
MR ratio (using the pessimistic ratio) so that R = AR is also advantageous regarding 
Johnson and shot noise providing a 10 — 20 dB boost in SNR (as explained above in section 

EC]). 

In addition the maximum MR ratio property allows to satisfy both requirements sought 
in CPP spin valves: small enough R A < 0.1 flfim 2 for high data rates, but large enough 
AR A > 5 mQ^im 2 for good SNR. (The latter figure would need to be adjusted to take into 
account the peculiarities of the non-local device but we use it as a reasonable rule-of-thumb). 

The main constraints for a large SNR will remain mag-noise, amplifier noise as well 
as the ability to scale down the non-local setup to par with the 20 — 30 nm track width 
relevantr^ 1 ^ for areal densities above 1 Tbit/inch 2 . Although we expect STT noise to be 
less of a nuisance, angular momentum transfer still occurs from injectors to detectors and a 
quantification of critical currents adapted to the geometry would be interesting. The angular 
response is thus an important aspect to explore in view of a potential use as a sensor and is 
discussed elsewhere 36 . 

Although the geometry studied in this paper might seem quite peculiar, there is actually 
a lot of flexibility in terms of design, provided the following recommendations are enforced: 

(i) confine spin with large interface resistances (when compared with the paramagnet 
spin resistance) wherever spin current can leak; 

(ii) in the smallest volume for the paramagnet connecting the ferromagnetic terminals 
(when compared to l^); 

(iii) to ensure infinite MR (optimistic) ratio, use two identical ferromagnetic electrodes 
as detectors within a symmetric arrangement with respect to injectors (the geometry of the 
ferromagnetic electrodes need not be completely identical: they just need to be identical to 
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within a few spin diffusion lengths lp of the contact to the paramagnet; this is an important 
practical remark since this allows to have different coercive fields for the various electrodes 
while keeping the properties stemming from symmetry); 

(iv) for further gain in the signal, use two antiparallel ferromagnetic electrodes. 

With these prescriptions non-local devices will perhaps be able to reach the realm of 
applications. 

V. CONCLUSION AND PROSPECTS. 

We have presented a variant of non-local spin valves based on a Van der Pauw cross 
setup: its basic characteristic is its reliance on four collinear ferromagnetic terminals. By 
playing on the numerous magnetization configurations of the four ferromagnets and on their 
symmetries several functionalities have been described: (i) an improved non-local spin valve 
(when used to read 1 bit) with an infinite GMR for the non-local resistance; (ii) ON-OFF 
switch effect where magnetization at the injectors can control the read-off at detectors; (iii) 
3-bit storage or sensing and offset-free 2-bit reading; (iv) direct Spin Hall measurement; (v) 
use as programmable magneto-logic gates. 

In addition to these functionalities the amplitude of the non-local resistance has the 
potential to be much increased in such a setup due to the conjunction of (i) spin confinement 
by tunnel barriers; (ii) the use of two electrodes as spin injectors instead of only one. We 
have also studied the separate or combined impact of both features by considering setups 
with two or three ferromagnets: the use of two injectors basically doubles the signal while the 
use of tunnel barriers hinder spin leaking from the Van der Pauw cross, thus increasing spin 
accumulation within it. The latter effect is already observed even when not all terminals are 
connected by tunnel barriers to the Van der Pauw cross. An important parameter is the spin 
resistance mismatch at each ferromagnetic-paramagnetic interface (at the four terminals), 
which is the ratio between the spin resistances of the ferromagnet and the paramagnet. The 
larger this parameter, the better the spin confinement in the Van der Pauw cross. 

Additional functionalities will be discussed elsewhere: non-collinear properties which are 
essential for sensor applications; as well as caloritronic ones. A related non-local device in 
the lateral geometry for which signals are much stronger will be also investigated^. 
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Appendix A: Computing the Spin Voltage. 

We derive in this Appendix the general expression of the spin voltage by relying on the 
one-dimensional drift-diffusion equations^i^. 



1. Basic equations. 

a. Spin accumulation vectors. 

We orient spin currents away from the origin O of the cross. 

In each electrode (Fl — FA) the spin accumulation and spin currents are then: 

z T 

A/i F1 (z) = Afx F1 (A) exp ( Al ) 

IsMz) = -Pnh - exp -if^ (A2) 

Kpi iFl 



so that: 



UA)=-W.-^- (A3) 



F2: without loss of generality we assume F2 antiparallel to Fl which allows to have 
symmetric equations (for the parallel case one should take opposite polarizations) 

W*) = -W« - ^f^ 1 (A4) 

where I Sy F2{B) is counted positive when flowing in the direction z' = —z; 
F3: 

WO = (A5) 

J*F3 

where I St Fs(B) is oriented in the direction x. 
FA: 

IsM D ) = 5 ( A6 ) 

where I s ,n^{D) is oriented in the direction x' = —x. 

In the arms (I — IV) we define spin accumulation vectors which are related to the spin 
accumulation and spin currents by: 
I: 



57 



Afn(z) = K : . (exp -i; exp ~ j , (A7) 
V In In / 

hA z ) = • ( ex P t-; - ex P -j- J (A8) 
Kjv \ In In/ 



where Ki is a constant two-component vector. 

II: following our convention of orienting away from the origin 



Afi H {z') = K„ . (exp i- ; exp -f - j , (A9) 
\ In In/ 

Is,n(z') = ^ • f exp ^; - exp -f] (A10) 



III: 



IV: 



fx X \ 

Afi in (x) = Km . exp — ; exp -— , (All) 
\ In In/ 

/ \ K m fx x \ /At n\ 

I s ,ni{x) = -— ■ exp — ; - exp -— (A12) 
Kn \ In In / 



fx' x'\ 
Afj, IV (x') = Kiv • exp — ; exp -— , (A13) 
V In In/ 

I s ,ivW) = ^ • fexp f; - exp -f) (A14) 
Kn V In In / 



b. Conditions at ferromagnetic - paramagnetic interfaces and at cross center. 

Interfaces. At interfaces we will neglect all spin flip but assume interface resistance is 
spin dependent. This implies that spin current at points A — D is continuous: 

I,, F1 (A) = I sJ (A) (A15) 

and so forth at points B — D. 

At points A — D the interface resistance induces a discontinuity in spin accumulations: 

AfiFi(A) - Am(A) = R cl I S {A) + R cl P cl I c (A16) 

Afi F2 {B) - A m (B) = R c2 I S {B) + R c2 P c2 h (A17) 
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Afi F3 (C) - Afi m {C) = R c3 I S {C) 
Afi F4 {D) - Afi IV (D) = R c4 I S {D) 

Cross center. At O, the spin accumulations are continuous: 

Am(0) = A/j n (B) = A mi {C) = Afi IV (0) 
and spin current continuity implies: 

o = i s>I (o) + i s>II (o) + LMO) + i.MO) 

so that for a = I — IV: 

K a . (1; 1) = A/i(0); 
£)K a .(l;-l) = 0. 

a 

2. Solution of the equations. 

a. Spin accumulation vectors K a 

We first solve for all four spin accumulation vectors in the arms (J — IV) 
From the equations ( 1A31 IA151 IA16I) at interface Fl - I one gets: 



A^(A) + (R cl + R F1 ) I S (A) = - (P F1 R F1 + P cl R cl ) I c 
After substitution in terms of spin accumulation vector Kj: 

Kj. ((1 + X X ) expZu (1 - X x ) exp-h) 

= — {Pf\Rf\ + PclRcl) Ic 



where: 



and 



R F i + R cl 

Ai — 



i? 



A' 



z - Ll 

In 
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Kn obeys a similar equation by changing the index 1 with index 2 since we have assumed 
that the magnetization of F2 is antiparallel to that of Fl. 



K„. ((1+X 2 ) exp/ 2 ; (1 - X 2 ) exp -l 2 ) 

- (P F2 R F 2 + Pc2Rc2) Ic (A28) 

while the equations for arms III — IV obtain by cancelling the charge current: 



Km • ((1 + X 3 ) exp l 3 ; (1 - X 3 ) exp -Z 3 ) = 0, (A29) 
K IV • ((1 + X A ) exp k; (1 - X A ) exp -Z 4 ) = 0. (A30) 

with X2/3/4 and Z2/3/4 defined as Xl and /1. 

Using the continuity of spin accumulation and spin current at O yields a second equation 
for each spin accumulation vector: 



K a . (1; l) = A/i(0), (a = I,.,IV) (A31) 

Solving for the K Q vectors in terms of Afx(0) leads to the following expressions: 

Afi(0) I (Xi - 1) exp -h 



while: 



where: 



K I = 



26 t \ (Xi + 1) exp h 



+ (lE^ELt^lM Ie ^ l (A32 ) 
K„ = 1 < — > 2 (A33) 



25 3 + ^ (X3 + I) exp/ 3 
Kiv = 3 < — > 4, (A35) 



± (Xj + 1) (X - 1) 
t> 4 = r exp/j ± exp — /j. (A36) 
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b. Spin currents at cross center. 



The relation between spin currents at O on each arm and the spin accumulation A/x(0) 
is derived using: 

LAO) = 4- K i ■ (!; - 1 ) ( A3r ) 

so that: 



_An(0) 

with a similar relation for arm II where: 



LAO) = —f^L - P em I c (A38) 

K eff,l 



5? 

Reff,i = Rn t^t (A39) 



P R 

PefU = (A40) 



5 

and: 



PRi = [P Fi R Fi + PciRa] /R N (A41) 
For arm III, one has simply: 

WO) = (A42) 

and a similar relation for arm IV. 

c. Spin accumulation at cross center. 

Plugging the expressions of the spin currents into the spin current continuity equation at 
O [Eq. (1A21I) ] allows determination of A/z(0): 



A//(0) = - (P eff>1 + P e//i2 ) 1— j- I c . (A43) 

Z-fi=l-4 - K e//,i 

Observe that the expression of A/j,(0) does not depend on the polarizations of the fer- 
romagnetic electrodes used as detectors (F 3 and F 4 ): it is fully determined by the injectors 
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only (Fl and F2). The previous expression has been derived by assuming Fl and F2 an- 
tiparallel, which reflects into the polarizations signs. In the more general case of arbitrary 
collinear orientation, one has: 

AMO) = - (P e //,1 " Pefffl) Reff U (A44) 

where positive polarizations are defined relative to an absolute axis (negative polarizations 
in the opposite direction). The largest spin accumulation in magnitude is achieved for 
antiparallel injector and collector electrodes. 

This can be plugged in the expressions for the spin current: 

(Peff,l ~ Peff,2) Reff 



ls,l(0) 



R 



effA 



eff,l 



(A45) 



with a similar expression for arm 77 while for arm 77/ (or IV with appropriate substitutions): 

(Peff,l - Pefffi) Reff 



IsJII 



(O) 



R 



eff,3 



Ic. 



(A46) 



d. Spin voltage. 



The voltage is measured as: 

V n i = - [/i F3 (+oo) - /i F4 (+oo)] /e. 



(A47) 



This voltage depends on the orientations of the four ferromagnetic terminals Fl — F4 so 
that in the most general case it can assume sixteen distinct values (2 4 ); but since the setup 
is invariant under a joint flipping of all the magnetizations this number is reduced to eight. 

We will need the chemical potential variation at interfaces: 



A*fs(C) - fim(C) = -R c3 P c3 7 a (C); 
H IV (D) - m (L>) = R c4 P c4 I S (D). 

Leaving out the factor e: 

-V n i = [^ 3 (+oo) - fi F3 {C)] + [fi F3 {C) - Hni(C)} 
+ [pi in (C) - pL in (0)} + [pu v (0) - fiiv(D)] 
+ \niv{D) - Hfa{D)\ + [fi F4 (D) - At F4 (+oo)] 
= P F3 An F3 (C)-R c3 P c3 I s (C) 
-P Fi A/i F4 ( J D) + R c4 P c4 I S (D) 
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so that: 



-V, 



nl 



PfzRfz + P&R, 



c3 



+ 



RF3 

PfaRfA + PciRci 
Rf4 

PfzRfz + P&Rc?, 

RF3 + Rc3 
PfaRfa + PcaRcA 



A» F3 (C) 
Afi F4 {D) 
A mi {C) 



RfA + Rc4 

The spin accumulations Ahjjj(C) and Afijy(D) are derived easily from the spin accu- 
mulation vectors Km and Kiv [Eq. (1A34HA35]) ]: 



Afi IH (C) = K m . (exp/ 3 ; exp -Z 3 ) 
A/i((9) 



*3 + 



(A48) 



and: 



Afi IV (D) = 3^4 



One thus gets a voltage drop: 



V nl = -A/i(0) 



PRa PRa 



5+ St 



-A/i((9) (P e// , 3 - PeZ/,4) 



(A49) 
(A50) 



e. Spin currents and spin accumulations at detector terminals. 



The spin current at point C is given by: 



so that: 



= -5- K m • (exp / 3 ; - exp -Z 3 ) (A51) 



= A /^(Q) = (Peff,l-Peff,2) Reffjc ^ 



Let us study |/ s (C)|as a function of X 3 while keeping all other parameters fixed. 
|/ s (C)|oc 1 



a6£ + S 3 

1 

oc 



2 { AT 3 [a cosh Z 3 + sinh Z 3 ] + [a sinh Z 3 + cosh Z 3 ] } 
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where: 

a = E Km > °- 

The spin current is therefore (in magnitude) a decreasing function of X 3 : in order to have 
large spin currents it is therefore better to have transparent junctions rather than large 
tunnel barriers. 

The spin accumulation at C behaves in an opposite manner: 

k 

oc 



aS 3 + 5 3 



X 3 

oc 



X 3 [a cosh l 3 + sinh / 3 ] + [a sinh l 3 + cosh l 3 ] 
which is an increasing function of X 3 . A large accumulation at the interface with detector 
F3 is favored by large spin impedance mismatch X 3 . 

The spin currents and accumulations at source and drain can be studied in the same 
manner with identical conclusions: large X leads to larger spin accumulations but smaller 
spin currents; with an opposite result for small X. 

3. Local resistance and non-local charge current. 

The resistance can also be measured locally along the charge current path from source 
to drain; if additionally we close the circuit between the two detector electrodes a charge 
current is generated. 

a. Local resistance at injectors. 

The voltage probes at source and drain are positioned at z\ and z 2 which we will assume 
to be very large (>> If% the spin relaxation lengths in the ferro magnets) . 

Vlocal = - \HFl(z\) - ^F2(Z 2 )] ■ (A53) 

(we have left out the factor e). 

We will need the dependence of the chemical potentials as well as the interface disconti- 
nuities: 
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VFi(z) = -p* F1 (1 - P 2 Fl ) I c z + d- P F1 Ap F1 ( Z ) 

PF2(z') = P* F2 (l- P| 2 ) Ic Z' + C 2 - P F2 Afl F2 ( Z ') 

where Ci_ 2 are constants and: 



p F1 (A) - fij(A) = -R cl (I c + P cl I s (A)) 
p F2 (B) - p n {B) = R c2 (J c - P c2 I s {B)) 



Since: 



- Vi oca i = [p F i{zi) - Pfi(A)] + [pfi(A) - m{A)\ 
+ [pi(A) - p n (B)} + [p n (B) - p F 2{B)\ 

+ \pL F2 (B) - (J,F2(Z 2 )] 



eventually: 



Viocai 



PfiRfi + -Pel -Rcl 

R F i + R c i 

PF2RF2 + Pc2Rc2 



RF2 + Rc2 



Apj(A) 
Ap n (B) 



Ro Ic 



where R ({P\~ 2 }) is an even function of polarizations {P\- 2 } = Pfi, P C i, 
Ro {P1-2} = P*Fi (! - P Fi) *i + Pn k 



i=l-2 



+R C i+ 



P Fi R ci RFi — P C iR C i — 2P Fi P ci R F iR C i\ 



Rfi + Rci 



(It can be checked that R > as it should be.) 
Since: 



Am{A) = Kj . (exp/ i; exp ~h) = A ^ X t 



St 



(with a similar expression for Afiu(B)) 



V local = Ap(0) [P em - P eff>2 ] - R I C 

= -Reff [Peff,l ~ Peff,2} 2 h ~ RqI C 
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The resistance variation when one flips terminal F2 is therefore: 

Vap - V P 

ARlnral =RaP ~ Rp= ; 



Hocal —J^AP 

= 4Pe//,l Peff,2Reff (A62) 
b. Local resistance at detectors and non-local charge current. 

The previous expressions will help us to derive the expression of the non-local charge 
current flowing through the detectors F3 and FA when they are in closed circuit instead of 
being in open circuit. 

Suppose one drives a current V c through terminals F3 and FA but not through Fl and 
F2 exchanging the roles of injectors and detectors (J c = 0); 

Vu = Reff [Pe//,3 - PeffAf I' e + ^0 {Ps-*} I' e (A63) 

We can define the local conductance at detectors: 

I' c = 934V U (A64) 

where 

g 3A = [Reff [P e//i3 - P effA ] 2 + Ro {P3-4}] _1 • (A65) 
We now switch on the charge current I c at source and drain Fl and F2; when V c = 0, 

V34 = Vnl — Rnl h 

= R eff (P e ff,l ~ Peff,l) ( P eff,3 ~ P eff,4) l c- (A66) 

When both I c and I' c are switched on, by superposition one gets the general expression of 
the current flowing through detector electrodes: 

I'c = 934 (V u - R nl I c ) . (A67) 
If we short F3 and FA, V 34 = so that the non-local charge current is: 

Inl = -934 Rnl h (A68) 

or: 

T R eff (Pem ~ P eff,2) ( P eff,3 ~ Peff,A> r ( \ 

nl = R VP P 12 , o rp 1 Ic ( M9 > 

Reff [Peff,3 ~ PeffM + K { P 3-4I 
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F1 




F 4t 



III 



z" 



Figure 24. Left: geometry used by Hershfield-Zhao and Fert-Lee for the bipolar spin switch calcu- 
lations. Right: geometry in our simplified one-dimensional treatment which includes spin leaking 
to the bottom. 

Appendix B: Revisiting the bipolar spin switch. 

The original theory by Johnson of the bipolar spin switch was refined by Hershfield and 
Zhao 48 who included spin relaxation in the ferromagnets, and also by Fert and Lee 4 ^ who 
studied the (adverse) impact of interface spin flips. If we discard the effect of spin flips, the 
non-local resistance variation was found by Hershfield and Zhao to be: 



where we have renamed parameters according to the definitions used in this paper: X 



and PR = (P C R C + PfRf) / Rn are defined at injector while Y and PR pertain to the 



detector. Fert and Lee have the same expression in the case PR = PR and X = Y (identical 
ferromagnets at injector and detector). 

What is noteworthy is that this expression is actually identical to that we found for 
the four ferromagnetic terminal cross geometry. In the limit of large X and Y , Ai? n ; 
might therefore be quite large as has been discussed repeatedly in this paper due to spin 
confinement. 

But a look at the geometry may cast doubt on this result (see Fig. [21]): there should 
be some spin leakage in the setup even if both injector and detector are in tunnel contacts 
with the central paramagnet, because spin is leaving along with charge on the bottom 



^Rnl 



AR n PR PR 



(Bl) 



[(X + 1) (Y + 1) exp / - (X - 1) (Y - 1) exp -/] 
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paramagnetic arm used as (charge) current drain. Hershfield and Zhao argued that spin 
leakage could be neglected in the bottom arm since in Johnson's original experiment the 
voltage pads are much larger than the spin relaxation length (the pads area was 0.01 mm 2 ^> 
l 2 N ) but it is difficult to see how this can be relevant for the bottom arm which is the current 
drain. The argument however applies to the spin leaking through the voltage probes. 

An expression which takes into account spin flow in the bottom arm is actually easy to 
derive in line with the calculations done with the cross geometry. 

Referring to Fig. [M] we consider a three arm geometry (arms / — III) with two ferro- 
magnetic terminals Fl and F2 with current flowing from Fl to III through /. We allow a 
finite size for arms / and II (L\ — 1\In and L 2 = I2 In) for a more general expression. We 
find: 



AR 



III 



2R eff P eff,l Peff,2 



(B2) 



where: 



R eff 



1 1 1 
+ -= — + 



R N R eff,l R eff,2\ 



(B3) 



The effective polarizations P e ff,i for terminals F1/F2 and spin resistances R e ff,i for arms I 
and 77 are defined as in the bulk of the paper (with Xi = X and X 2 = Y): 



p— 1 _ p-i 

U eff - N 



R N <5j #2 

R' * + + 



(B4) 



„ N 6t 5 2 + 

We have allowed for different cross sections An and A' N along x'Ox and Oz' so that the 
spin resistance along x'Ox is Rn = p* n In/An but along Oz' is 

r 'n — P*nIn/A' n . 

The total spin resistance is easy to interpret: it corresponds to parallel addition of the spin 
resistance of the three arms with spin resistance R' N for arm 777 and spin resistances R e ff,i 
and R e ff, 2 for arms I and II respectively. 

Suppose we neglect arm III spin resistance so that: 

1 



R„tt 



£ff R eff,l + R eff,2 

One then recovers Fert-Lee and Hershfield-Zhao expression [Eq. (jBl[) ] when l\ = l 2 = 1/2 
which therefore does imply a neglect of spin leakage to the bottom terminal. 



(B5) 
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When is it allowed to do so? The condition is: 



R' N » (R e ff t i, R e ff,2) (B6) 
The effective resistance varies between Rn and Rp + R c ; the condition then reduces to: 

R N S> (Rfi + -Rci, Rf2 + R C 2) Rn) ■ (B7) 

which is favored by small cross-section <C • Note that the condition R^ <C R' N is 
however not enough to recover Eq. ( IB II) . 

So more generally when the condition in Eq. flB7j) is not met, one gets (when l\ = I2 = 

1/2): 

Ai^ = 8R n PR PR 

x {2 (X + 1) (Y + 1) exp I - 2 (X - 1) (Y - 1) exp -/ 

+ [( X + 1)(Y+ 1) exp I + (X - 1) (Y - 1) exp -/ + 2 (XF - l)]}- 1 (B8) 

which is significantly different from Eq. (jBlft . Let's define the spin mismatch Z = R' n /Rn- 
At small distance where the signal will be largest: 

XY 



AR 



ni 



X- 1 + Y- 1 + Z- 1 ' 

The scale will be set by the smaller of the spin impedance mismatches. If anyone of them 
is in the transparent regime, one will recover the scaling characteristic of open systems: 
^Rni < -Rat as it should be. 

We have neglected in the above spin leakage in the voltage probes but it is easy to 
generalize the present considerations to include it, should it become experimentally relevant. 
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